Biomarker Panel for Diagnosis and Prediction of Graft Rejection

ABSTRACT

Methods are provided for monitoring a subject having a graft for an acute rejection (AR) response, e.g., to predict, to diagnose, and/or to characterize an AR response. In practicing the subject methods, the expression level of at least one gene in a sample from the subject, e.g., a blood or biopsy sample, is evaluated, e.g., at the nucleic acid and/or protein level, to monitor the subject. Also provided are compositions, systems, kits and computer program products that find use in practicing the subject methods.

GOVERNMENT RIGHTS

This invention was made with Government support under contract (R01 AI 61739-01) awarded by the National Institutes of Health. The Government has certain rights in this invention.

BACKGROUND

Transplantation of a graft organ or tissue from a donor to a host patient is a feature of certain medical procedures and treatment protocols. Despite efforts to avoid graft rejection through host-donor tissue type matching, in transplantation procedures where a donor organ is introduced into a host, immunosuppressive therapy is generally required to the maintain viability of the donor organ in the host. However, despite the wide use of immunosuppressive therapy, organ transplant rejection can occur.

Acute graft rejection (AR) of allograft tissue is a complex immune response that involves T-cell recognition of alloantigen in the allograft, co-stimulatory signals, elaboration of effectors molecules by activated T cells, and an inflammatory response within the graft. Activation and recruitment of circulating leukocytes to the allograft is a central feature of this process.

Early detection of AR is one of the major clinical concerns in the care of transplant recipients, including kidney transplant recipients. Detection of AR before the onset of renal dysfunction allows successful treatment of this condition with aggressive immunosuppression. It is equally important to reduce immunosuppression in patients who do not have AR to minimize drug toxicity.

Accordingly, techniques for monitoring for an AR response in a transplant recipient, including predicting, diagnosing and characterizing AR, are of interest in the field. The present invention meets these and other needs.

SUMMARY OF THE INVENTION

Methods are provided for monitoring a subject having a graft for an acute rejection (AR) response, e.g., to predict, to diagnose, and/or to characterize an AR response, including graft survival in a subject. In practicing the subject methods, the expression level of at least one gene in a sample from the subject, e.g., a blood or biopsy sample, is evaluated, e.g., at the nucleic acid and/or protein level, to monitor the subject. Also provided are compositions, systems, kits and computer program products that find use in practicing the subject methods. The methods and compositions find use in a variety of applications.

Aspects of the subject invention include methods of monitoring a subject who has received a graft (e.g., an allograft) for an acute rejection (AR) response including: (a) evaluating the expression level of at least one gene in a sample from the subject to obtain a gene expression result, wherein the at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 (also called NAMPT) and PSEN1; and (b) employing the gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in the subject, thereby monitoring the subject for an AR response. In certain embodiments, the expression level of all of the genes from the above 5-gene group is assessed.

In certain embodiments, the evaluating step further includes evaluating the expression level of additional genes selected from one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP. In certain embodiments, the expression level of all of the 10 genes listed above is assessed. In certain embodiments, the evaluating step comprises assaying the sample for an expression product of the at least one gene. In certain embodiments, the expression product is a nucleic acid transcript whereas in other embodiments, the expression product is a protein. In certain embodiments, the assessing step comprises an RT-PCR assay. In certain embodiments, the RT-PCR assay is a quantitative RT-PCR assay. In certain embodiments, the sample is a blood sample. In certain embodiments, the allograft is a kidney allograft. In certain embodiments, the gene expression result is employed to predict the occurrence of an AR response between 6 and 3 months in advance. In certain embodiments, the gene expression result is employed to characterize an AR response as a steroid resistant AR response.

Aspects of the invention include methods for determining the identity of a tissue that is being rejected in an AR response, including evaluating the expression level of at least one gene selected from one or more of NKTR, MAPK9, DUSP1, PBEF1 (also called NAMPT), PSEN1, CFLAR, RNF-130, IFNGR1, ITGAX and RYBP in a sample from the subject to obtain a gene expression result and employing the gene expression result to determine the identity of a tissue that is being rejected. The test is applicable not just transplant rejection of numerous tissues, including kidney, heart, liver, lung and intestinal transplants.

Aspects of the subject invention include methods for determining an immunosuppressive regimen for a subject who has received a graft including: (a) evaluating the expression level of at least one gene in a sample from the subject to obtain a gene expression result, wherein the at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 (also called NAMPT) and PSEN1; and (b) employing the gene expression result to determine an immunosuppressive regimen for the subject. In certain embodiments, the expression level of all of the genes from the above 5-gene group is assessed. In certain embodiments, the evaluating step further includes evaluating the expression level of additional genes selected from one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP. In certain embodiments, the expression level of all of the 10 genes listed above is assessed.

Aspects of the subject invention include systems for monitoring a subject who has received an allograft for an acute rejection (AR) response, the system including: (a) a gene expression evaluation element configured for evaluating the expression level of at least one gene in a sample from the subject to obtain a gene expression result, wherein the at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 (also called NAMPT) and PSEN1; and (b) a phenotype determination element for employing the gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in the subject, thereby monitoring the subject for an AR response. In certain embodiments, the expression level of all of the genes from the above 5-gene group is assessed. In certain embodiments, the gene expression evaluation element comprises at least one reagent for assaying a sample for an expression product of the at least one gene. In certain embodiments, the expression product of the at least one gene is a nucleic acid transcript whereas in other embodiments, the expression product is a protein. In certain embodiments, the gene expression evaluation element comprises PCR primers specific for the at least one gene. In certain embodiments, the gene expression evaluation element is further configured for evaluating the expression level of additional genes selected from one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP. In certain embodiments, the expression level of all of the 10 genes listed above is assessed. In certain embodiments, the phenotype determination element comprises a reference expression value for the at least one gene.

Aspects of the subject invention include kits for monitoring a subject who has received an allograft for an acute rejection (AR) response, including: (a) a gene expression evaluation element for evaluating expression of at least one gene in a sample to obtain a gene expression result, wherein the at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 (also called NAMPT) and PSEN1; and (b) a phenotype determination element for employing the gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in the subject, thereby monitoring the subject for an AR response. In certain embodiments, the expression level of all of the genes from the above 5-gene group is assessed. In certain embodiments, the gene expression evaluation element is further configured for evaluating the expression level of additional genes selected from one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP. In certain embodiments, the expression level of all of the 10 genes listed above is assessed. In certain embodiments, the gene expression evaluation element includes PCR primer pairs specific for one or more of the genes noted above (e.g., for use in an RT-PCR assay for gene expression). In certain embodiments, the gene expression evaluation element includes PCR primer pairs specific for all of the 5 genes in the 5 gene group, whereas in other embodiments, the gene expression evaluation element includes PCR primer pairs specific for all 10 genes listed above.

Aspects of the subject invention include computer program products for monitoring a subject who has received an allograft for an acute rejection (AR) response, wherein the computer program product, when loaded onto a computer, is configured to employ a gene expression result from a sample derived from the subject to determine an AR monitoring result and provide the AR monitoring result to a user in a user-readable format, wherein the AR monitoring result is selected from predicting the onset of an AR response, diagnosing an AR response, and/or characterizing an AR response in the subject, wherein the gene expression result includes expression data for one or more genes selected from: NKTR, MAPK9, DUSP1, PBEF1 and PSEN1. In certain embodiments, the expression result includes expression data for all of the genes from the above 5-gene group. In certain embodiments, the expression result includes expression data for additional genes selected from one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP. In certain embodiments, the expression result includes expression data for all of the 10 genes listed above.

Definitions

For convenience, certain terms employed in the specification, examples, and appended claims are collected here.

“Acute rejection or AR” is the rejection by the immune system of a tissue transplant recipient when the transplanted tissue is immunologically foreign. Acute rejection is characterized by infiltration of the transplanted tissue by immune cells of the recipient, which carry out their effector function and destroy the transplanted tissue. The onset of acute rejection is rapid and generally occurs in humans within a few weeks after transplant surgery. Generally, acute rejection can be inhibited or suppressed with immunosuppressive drugs such as rapamycin, cyclosporin A, anti-CD40L monoclonal antibody and the like.

“Chronic transplant rejection or CR” generally occurs in humans within several months to years after engraftment, even in the presence of successful immunosuppression of acute rejection. Fibrosis is a common factor in chronic rejection of all types of organ transplants. Chronic rejection can typically be described by a range of specific disorders that are characteristic of the particular organ. For example, in lung transplants, such disorders include fibroproliferative destruction of the airway (bronchiolitis obliterans); in heart transplants or transplants of cardiac tissue, such as valve replacements, such disorders include fibrotic atherosclerosis; in kidney transplants, such disorders include, obstructive nephropathy, nephrosclerorsis, tubulointerstitial nephropathy; and in liver transplants, such disorders include disappearing bile duct syndrome. Chronic rejection can also be characterized by ischemic insult, denervation of the transplanted tissue, hyperlipidemia and hypertension associated with immunosuppressive drugs.

The term “transplant rejection” encompasses both acute and chronic transplant rejection.

The term “stringent assay conditions” as used herein refers to conditions that are compatible to produce binding pairs of nucleic acids, e.g., surface bound and solution phase nucleic acids, of sufficient complementarity to provide for the desired level of specificity in the assay while being less compatible to the formation of binding pairs between binding members of insufficient complementarity to provide for the desired specificity. Stringent assay conditions are the summation or combination (totality) of both hybridization and wash conditions.

“Stringent hybridization conditions” and “stringent hybridization wash conditions” in the context of nucleic acid hybridization (e.g., as in array, Southern or Northern hybridizations) are sequence dependent, and are different under different experimental parameters. Stringent hybridization conditions that can be used to identify nucleic acids within the scope of the invention can include, e.g., hybridization in a buffer comprising 50% formamide, 5×SSC, and 1% SDS at 42° C., or hybridization in a buffer comprising 5×SSC and 1% SDS at 65° C., both with a wash of 0.2×SSC and 0.1% SDS at 65° C. Exemplary stringent hybridization conditions can also include hybridization in a buffer of 40% formamide, 1 M NaCl, and 1% SDS at 37° C., and a wash in 1×SSC at 45° C. Alternatively, hybridization to filter-bound DNA in 0.5 M NaHPO₄, 7% sodium dodecyl sulfate (SDS), 1 mM EDTA at 65° C., and washing in 0.1×SSC/0.1% SDS at 68° C. can be employed. Yet additional stringent hybridization conditions include hybridization at 60° C. or higher and 3×SSC (450 mM sodium chloride/45 mM sodium citrate) or incubation at 42° C. in a solution containing 30% formamide, 1M NaCl, 0.5% sodium sarcosine, 50 mM MES, pH 6.5. Those of ordinary skill will readily recognize that alternative but comparable hybridization and wash conditions can be utilized to provide conditions of similar stringency.

In certain embodiments, the stringency of the wash conditions that set forth the conditions which determine whether a nucleic acid is specifically hybridized to a surface bound nucleic acid. Wash conditions used to identify nucleic acids may include, e.g.: a salt concentration of about 0.02 molar at pH 7 and a temperature of at least about 50° C. or about 55° C. to about 60° C.; or, a salt concentration of about 0.15 M NaCl at 72° C. for about 15 minutes; or, a salt concentration of about 0.2×SSC at a temperature of at least about 50° C. or about 55° C. to about 60° C. for about 15 to about 20 minutes; or, the hybridization complex is washed twice with a solution with a salt concentration of about 2×SSC containing 0.1% SDS at room temperature for 15 minutes and then washed twice by 0.1×SSC containing 0.1% SDS at 68° C. for 15 minutes; or, equivalent conditions. Stringent conditions for washing can also be, e.g., 0.2×SSC/0.1% SDS at 42° C.

A specific example of stringent assay conditions is rotating hybridization at 65° C. in a salt based hybridization buffer with a total monovalent cation concentration of 1.5 M (e.g., as described in U.S. patent application Ser. No. 09/655,482 filed on Sep. 5, 2000, the disclosure of which is herein incorporated by reference) followed by washes of 0.5×SSC and 0.1×SSC at room temperature.

Stringent assay conditions are hybridization conditions that are at least as stringent as the above representative conditions, where a given set of conditions are considered to be at least as stringent if substantially no additional binding complexes that lack sufficient complementarity to provide for the desired specificity are produced in the given set of conditions as compared to the above specific conditions, where by “substantially no more” is meant less than about 5-fold more, typically less than about 3-fold more. Other stringent hybridization conditions are known in the art and may also be employed, as appropriate.

As used herein, the term “gene” or “recombinant gene” refers to a nucleic acid comprising an open reading frame encoding a polypeptide, including exon and (optionally) intron sequences. The term “intron” refers to a DNA sequence present in a given gene that is not translated into protein and is generally found between exons in a DNA molecule. In addition, a gene may optionally include its natural promoter (i.e., the promoter with which the exons and introns of the gene are operably linked in a non-recombinant cell, i.e., a naturally occurring cell), and associated regulatory sequences, and may or may not have sequences upstream of the AUG start site, and may or may not include untranslated leader sequences, signal sequences, downstream untranslated sequences, transcriptional start and stop sequences, polyadenylation signals, translational start and stop sequences, ribosome binding sites, and the like.

A “protein coding sequence” or a sequence that “encodes” a particular polypeptide or peptide, is a nucleic acid sequence that is transcribed (in the case of DNA) and is translated (in the case of mRNA) into a polypeptide in vitro or in vivo when placed under the control of appropriate regulatory sequences. The boundaries of the coding sequence are determined by a start codon at the 5′ (amino) terminus and a translation stop codon at the 3′ (carboxy) terminus. A coding sequence can include, but is not limited to, cDNA from viral, procaryotic or eukaryotic mRNA, genomic DNA sequences from viral, procaryotic or eukaryotic DNA, and even synthetic DNA sequences. A transcription termination sequence may be located 3′ to the coding sequence.

The terms “reference” and “control” are used interchangeably to refer to a known value or set of known values against which an observed value may be compared. As used herein, known means that the value represents an understood parameter, e.g., a level of expression of a marker gene in a graft survival or loss phenotype.

The term “nucleic acid” includes DNA, RNA (double-stranded or single stranded), analogs (e.g., PNA or LNA molecules) and derivatives thereof. The terms “ribonucleic acid” and “RNA” as used herein mean a polymer composed of ribonucleotides. The terms “deoxyribonucleic acid” and “DNA” as used herein mean a polymer composed of deoxyribonucleotides. The term “mRNA” means messenger RNA. An “oligonucleotide” generally refers to a nucleotide multimer of about 10 to 100 nucleotides in length, while a “polynucleotide” includes a nucleotide multimer having any number of nucleotides.

The terms “protein” and “polypeptide” used in this application are interchangeable. “Polypeptide” refers to a polymer of amino acids (amino acid sequence) and does not refer to a specific length of the molecule. Thus peptides and oligopeptides are included within the definition of polypeptide. This term does also refer to or include post-translational modifications of the polypeptide, for example, glycosylations, acetylations, phosphorylation and the like. Included within the definition are, for example, polypeptides containing one or more analogs of an amino acid, polypeptides with substituted linkages, as well as other modifications known in the art, both naturally occurring and non-naturally occurring.

The term “assessing” and “evaluating” are used interchangeably to refer to any form of measurement, and includes determining if an element is present or not. The terms “determining,” “measuring,” “assessing,” and “assaying” are used interchangeably and include both quantitative and qualitative determinations. Assessing may be relative or absolute. “Assessing the presence of” includes determining the amount of something present, as well as determining whether it is present or absent.

Certain abbreviations employed in this application include the following:

AR: Acute Rejection; AZA: Azathioprine; BKV: BK Virus; CMV: Cytomegalovirus; CSA: cyclosporine A; DSA: Donor Specific Antibodies; EBV: Epstein-Barr virus; FDR: false discovery rate; FK: FK506; HCT: Hematocrit; LRD: living related donors; MMF: Mycophenolate Mofetil; PBL: Peripheral Blood Leukocytes; PAM: Prediction Analysis of Microarrays; PRA: Panel Reactive Antibody; Q-PCR: quantitative real time polymerase chain reaction; ROC: Receiver Operating Characteristic; SAM: Significance Analysis of Microarrays; STA: stable; SNS: patients from the NIH-CCTPT U01 funded randomized multicenter pediatric renal transplant steroid-based vs. steroid-free transplantation study (SNS01); WBC: White blood cell.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1. Summary of Study Design. The present study was divided into 2 specific parts: In part 1 of the study: we conducted cross sectional microarray analysis of peripheral blood samples on 3 different microarray platforms—the Affymetrix, the Agilent and the cDNA arrays. Array data generated from the 3 platforms were compared by mapping the transcripts to Human Gene Organization (HUGO) gene names. Appropriate normalization procedures were applied to each dataset followed by SAM analysis for identification of genes regulated significantly in AR, using a common significance threshold (FDR<10%). 525 genes were common among all differentially expressed AR specific genes across the 3 platforms. A 32 gene-set was selected for initial verification on 34 samples (Verification Set) chosen from the samples used on the microarrays. Of 32 genes tested, the 10 most significant genes (p<0.03) were further validated on 42 independent samples (Validation Set1). A regression model was generated using the expression of 5 most significant genes from Validation Set 1. This model was then applied to the second independent set of 64 samples (Validation Set 2) for true AR classification. In Part 2 of the study, we performed longitudinal of 5-gene Q-PCR analysis of paired blood samples collected before (preAR; n=27) and after AR (postAR; n=30), in which the regression model was applied to 57 serial AR samples.

FIG. 2. Verification and Validation of Gene Expression for the 5 Gene-Set. Clustering of arrays from peripheral blood samples across the common 525 gene-set shows clustering of samples by diagnosis phenotype. Clustering of 59 samples run on Affymetrix arrays is shown here (FIG. 2A); clustering of samples by phenotype is seen across the 525 overlapping genes, differentially expressed on all 3 array platforms (see FIG. 5), irrespective of sample collection, processing or multiple clinical confounders. The samples cluster by phenotype (AR or STA) across the 525 genes. Shown here, are representative clustering data for the samples run only on the Affymetrix array. Summarized expression of the selected 10 gene-set is shown as measured on one of the array platforms (Affymetrix) (FIG. 2B). Fold changes and q score from SAM analysis of microarray data from Affymetrix in whole blood samples are shown. The expression of the same 10 genes is shown for the Verification Set of Q-PCR samples (FIG. 2C). The fold change in Q-PCR was the mean differences between AR and STA samples; the significance for expression differences by Q-PCR were determined by T test, where p<0.05 was taken to be significant. In the independent validation set of blood samples (Validation Set 1), the expression for the most significant 5 genes is shown. The new blood samples, not used in the array analyses, again segregate by phenotype (AR or STA) across the 5 genes (FIG. 2D). The fold changes and p values of the 5 genes in the samples used in Validation Set 1 are shown in FIG. 2E.

FIG. 3. Receiver Operating Characteristic (ROC) Curve. Using logistic regression model analysis, a prediction model was generated from 5 gene Q-PCR expression data of independent samples (Validation Set1). The specificity and sensitivity of the model is shown by the ROC score=95.2%. The regression model was applied to an independent sample set (Validation set 2), where the expression of 5 genes clearly separates AR from STA (unsupervised clustering), and the differentially expressed 5 genes are significantly different in AR when comparing to STA (p<0.02).

FIG. 4. Prediction of Acute Rejection by Blood Sampling for the 5 Gene-Set Prior to Clinical Graft Dysfunction. The regression model was applied to samples sequentially collected from patients in the STA sample set (FIG. 4A), at various times in the first 2 years post-transplantation, to evaluate the probability of prediction as AR at different times after transplantation. The STA samples were grouped to samples collected within 0 to 3 months, within 3 to 6 months, within 6 to 12 months and within 12 to 24 months post transplantation. Each of the blood samples had a paired biopsy to confirm absence of AR in the graft at the time of blood sample collection. The mean and standard error of AR probability was calculated for all STA samples using the coordinated expression of the 5 gene-set at different times after kidney transplantation. As seen, the likelihood of any of the STA samples being predicted as AR, at any time post-transplantation, was less than 50%. There were 2 misclassified STA samples (S99 with an AR probability score of 58% and S88 with an AR prediction score of 72%; see FIG. 9). Both patients with misclassified STA samples have gone on to develop biopsy proven AR within 6 months of the blood sampling. In the right panel (FIG. 4B), the mean and standard error of AR probability is shown for blood samples collected at the time of biopsy proven AR, as well as for samples collected (as part of protocol sample collection) within 1-3 months prior and 3-6 months prior to the AR sample (preAR). As seen, the probability of classification of the blood sample as AR (generally 75-100%), is in strong agreement with the biopsy diagnosis of AR. After intensification of immunosuppression at AR diagnosis (steroid pulsing or antibody), the AR prediction probability falls, demonstrating that the AR prediction probability of the 5 gene-set signature is reduced by augmenting immunosuppression. The AR prediction probability for all samples in the AR set, is significantly higher at all preAR, AR and postAR samples, than any of the STA samples (p<0.001), suggesting that patients who develop AR may have a higher immune activation threshold than patients who never develop AR after transplantation.

FIG. 5 provides Venn diagrams showing the significant genes for AR within each discovery platform and the overlapping significant genes in bold font, as ascertained by SAM (q<10%) across sample preparations (A) or across different microarray platforms (B).The overlapping 525 genes (double underline) are statistically significantly (Hypergeometric test; p<0.001) regulated in AR, compared to STA samples and are common to AR, irrespective of hybridization to different array platforms are the use of different methods of peripheral blood collection.

FIG. 6A shows a list of genes and ABI probe ID for genes assayed using Q-PCR. FIG. 6B shows normalized expression values for 10-genes by Q-PCR across 34 blood samples in the verification set (17 AR, 17STA).

FIG. 7 is a summary of gene expression values on array and Q-PCR across the 34 samples in the verification set for 3 genes indicated in the literature as differentially expressed in AR.

FIG. 8 shows normalized expression values for 10 genes across 42 independent blood samples (21 AR and 21 STA) in validation set 1 by Q-PCR.

FIG. 9 shows normalized expression values for 5 genes across 64 independent blood samples (32 AR and 32 STA) in validation set 2 by Q-PCR.

FIG. 10 compares different classification models for AR across the 10 gene set (top panel) and 5 gene set (bottom panel).

FIG. 11 shows the 18 demographic and clinical Confounders examined for all Q-PCR experiments.

FIG. 12 shows results of Pearson correlation and t-test (as appropriate) to determine any influence of the 18 parameters in FIG. 11 as potential confounders for expression of the 5-genes set.

FIG. 13 is a graphical representation of the results in FIG. 12.

FIG. 14 shows the correlation between biopsy Banff grade, Peritubular C4d positivity and humoral vs. cellular AR and the AR prediction probability.

FIG. 15. Panel A shows prediction probability of AR in paired samples collected within 6 months before and after the AR episode. Panel B shows p values comparing the AR prediction score at AR vs. preAR and postAR scores. Panel C is a graphical representation of the data in panel A. Panel D is a graph showing AR probability in STA samples. Panel E shows the data used in the graph shown in panel D.

FIG. 16 shows biological relevance analysis for the overlapping 525 genes across the 3 array platforms using Ingenuity System (IPA; 2008).

FIG. 17 is a graphical representation of the information shown in FIG. 16.

FIG. 18 shows a table cross-referencing the samples and assays performed.

DESCRIPTION OF THE SPECIFIC EMBODIMENTS

Methods are provided for monitoring a subject having a graft for an acute rejection (AR) response, e.g., to predict, to diagnose, and/or to characterize an AR response. In practicing the subject methods, the expression level of at least one gene in a sample from the subject, e.g., a blood or biopsy sample, is evaluated, e.g., at the nucleic acid and/or protein level, to monitor the subject. Also provided are compositions, systems, kits and computer program products that find use in practicing the subject methods. The methods and compositions find use in a variety of applications.

Before the present invention is described in greater detail, it is to be understood that this invention is not limited to particular embodiments described, as such may, of course, vary. It is also to be understood that the terminology used herein is for the purpose of describing particular embodiments only, and is not intended to be limiting, since the scope of the present invention will be limited only by the appended claims.

Where a range of values is provided, it is understood that each intervening value, to the tenth of the unit of the lower limit unless the context clearly dictates otherwise, between the upper and lower limit of that range and any other stated or intervening value in that stated range, is encompassed within the invention. The upper and lower limits of these smaller ranges may independently be included in the smaller ranges and are also encompassed within the invention, subject to any specifically excluded limit in the stated range. Where the stated range includes one or both of the limits, ranges excluding either or both of those included limits are also included in the invention.

Certain ranges are presented herein with numerical values being preceded by the term “about.” The term “about” is used herein to provide literal support for the exact number that it precedes, as well as a number that is near to or approximately the number that the term precedes. In determining whether a number is near to or approximately a specifically recited number, the near or approximating unrecited number may be a number which, in the context in which it is presented, provides the substantial equivalent of the specifically recited number.

Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although any methods and materials similar or equivalent to those described herein can also be used in the practice or testing of the present invention, representative illustrative methods and materials are now described.

All publications and patents cited in this specification are herein incorporated by reference as if each individual publication or patent were specifically and individually indicated to be incorporated by reference and are incorporated herein by reference to disclose and describe the methods and/or materials in connection with which the publications are cited. The citation of any publication is for its disclosure prior to the filing date and should not be construed as an admission that the present invention is not entitled to antedate such publication by virtue of prior invention. Further, the dates of publication provided may be different from the actual publication dates which may need to be independently confirmed.

It is noted that, as used herein and in the appended claims, the singular forms “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise. It is further noted that the claims may be drafted to exclude any optional element. As such, this statement is intended to serve as antecedent basis for use of such exclusive terminology as “solely,” “only” and the like in connection with the recitation of claim elements, or use of a “negative” limitation.

As will be apparent to those of skill in the art upon reading this disclosure, each of the individual embodiments described and illustrated herein has discrete components and features which may be readily separated from or combined with the features of any of the other several embodiments without departing from the scope or spirit of the present invention. Any recited method can be carried out in the order of events recited or in any other order which is logically possible.

As summarized above, the subject invention is directed to methods of monitoring for an AR response in a subject, as well as reagents and kits for use in practicing the subject methods. In further describing the invention, the subject methods are described first, followed by a review of the reagents and kits for use in practicing the subject methods.

Methods of Monitoring for an Acute Rejection Response

As reviewed above, the subject invention provides methods for monitoring a subject who has received a graft for an AR response. In certain embodiments, a subject is monitored to predict if an AR response will occur, where the method can predict a AR from between 6 and 3 months prior to the occurrence of an AR response. In certain embodiments, a subject is monitored to determine whether there is an ongoing AR response. In certain embodiments, a subject is monitored to characterize an AR response that is or has occurred, e.g., to determine the severity and/or class of a previous AR response, e.g., whether an AR response is/was a steroid resistant AR response.

In certain embodiments, the subject invention includes methods for determining the identity of a tissue that is being rejected in an AR response. Aspects of the subject invention further include methods for determining an immunosuppressive regimen for a subject who has received a graft, e.g., an allograft.

As such, certain embodiments of the invention provide methods of evaluating, e.g., in terms of predicting, graft survival in a subject comprising a graft. For example, the subject methods can be used to determine whether a subject will mount an AR response to a graft that, in the absence of effective immunosuppressive intervention, will result in loss of the graft tissue. As such, the subject invention provides methods of evaluating whether a graft in a transplant patient or subject will survive or be lost. In certain embodiments, the methods may be viewed as methods of determining whether a transplant subject has a graft survival phenotype, i.e., a phenotype in which the graft will survive. A graft survival phenotype is a phenotype characterized by the presence of long-term graft survival. By “long-term” graft survival is meant graft survival for at least about 5 years beyond current sampling, despite the occurrence of one or more prior episodes of AR. In certain embodiments, graft survival is determined for patients in which at least one episode of acute rejection (AR) has occurred. As such, these embodiments are methods of determining or predicting graft survival following AR. Graft survival is determined or predicted in certain embodiments in the context of transplant therapy, e.g., immunosuppressive therapy, where immunosuppressive therapies are known in the art. In yet other embodiments, methods of determining the class and/or severity of acute rejection (and not just the presence thereof) are provided.

As in known in the transplantation field, a graft organ, tissue or cell(s) may be allogeneic or xenogeneic, such that the grafts may be allografts or xenografts. Organs and tissues of interest include, but are not limited to: skin, heart, kidney, liver, bone marrow, and other organs.

In practicing the subject methods, a subject or patient sample, e.g., cells or collections thereof, e.g., tissues, is assayed to monitor for an AR response in the host. Accordingly, the first step of the subject methods is to obtain a suitable sample from the subject or patient of interest, i.e., a patient having at least one graft, e.g., allograft.

The sample is derived from any initial suitable source, where sample sources of interest include, but are not limited to, many different physiological sources, e.g., cerebro-spinal fluid (CSF), urine, saliva, tears, tissue derived samples, e.g., homogenates (such as biopsy samples of the transplanted tissue or organ (including, but not limited to kidney, heart, lung biopsies)), and blood or derivatives thereof.

In certain embodiments, a suitable initial source for the patient sample is a blood sample. As such, the sample employed in the subject assays of these embodiments is generally a blood-derived sample. The blood-derived sample may be derived from whole blood or a fraction thereof, e.g., serum, plasma, cellular fraction, etc., where in certain embodiments the sample is derived from blood cells harvested from whole blood. Of particular interest as a sample source are peripheral blood mononuclear cells/lymphocytes (PBMCs/PBLs). Any convenient protocol for obtaining such samples may be employed, where suitable protocols are well known in the art (e.g., density gradient fractionation of a whole blood sample) and a representative protocol is reported in the Experimental Section, below.

In practicing the subject methods, the sample is assayed to obtain an expression level evaluation, e.g., expression profile, for one or more genes selected from Table 1, where the term expression profile is used broadly to include a genomic expression profile, e.g., an expression profile of nucleic acid transcripts, e.g., mRNAs, of the one or more genes of interest, or a proteomic expression profile, e.g., an expression profile of one or more different proteins, where the proteins/polypeptides are expression products of the one or more genes of interest. As such, in certain embodiments the expression level of only one gene in Table 1 is evaluated. In yet other embodiments, the expression level of two or more genes from Table 1 is evaluated, e.g., 3, 4 or all 5 genes in Table 1. In certain embodiments, the expression level of one or more additional gene other than those listed in Table 1 is also evaluated, where in certain embodiments, the additional gene(s) is selected from those listed in Table 2. It is noted here that an expression profile that includes an evaluation of the expression level of any combination of genes in Tables 1 and 2 finds use in monitoring a subject having or being predisposed to develop an AR response in the engrafted tissue, including evaluating the expression of all 10 genes listed in Tables 1 and 2.

TABLE 1 5 genes whose expression level can be used to monitor AR in a subject having an allograft. NCBI Gene Symbol Gene Name GeneID No. DUSP1 dual specific kinase 1 1843 MAPK9 mitogen-activated protein kinase 5601 NKTR natural killer-tumor recognition sequence 4820 PBEF1 pre-B cell colony enhancing factor 10135 (or NAMPT) PSEN1 presenilin 1 5663

TABLE 2 5 additional genes whose expression level can be used to monitor AR in a subject having an allograft. Gene NCBI Symbol Gene Name GeneID No. CFLAR caspase 8 and FADD-like apoptosis regulator 8837 IFNGR1 interferon gamma receptor 1 3459 ITGAX Integrin alpha X 3687 RNF130 ring finger protein 130 55819 RYBP ring 1 and YY1 binding protein 23429

In the broadest sense, the expression evaluation may be qualitative or quantitative. As such, where detection is qualitative, the methods provide a reading or evaluation, e.g., assessment, of whether or not the target analyte, e.g., nucleic acid or other expression product (e.g., protein), is present in the sample being assayed. In yet other embodiments, the methods provide a quantitative detection of whether the target analyte is present in the sample being assayed, i.e., an evaluation or assessment of the actual amount or relative abundance of the target analyte, e.g., nucleic acid in the sample being assayed. In such embodiments, the quantitative detection may be absolute or, if the method is a method of detecting two or more different analytes, e.g., target nucleic acids, in a sample, relative. As such, the term “quantifying” when used in the context of quantifying a target analyte, e.g., nucleic acid(s), in a sample can refer to absolute or to relative quantification. Absolute quantification may be accomplished by inclusion of known concentration(s) of one or more control analytes and referencing the detected level of the target analyte with the known control analytes (e.g., through generation of a standard curve). Alternatively, relative quantification can be accomplished by comparison of detected levels or amounts between two or more different target analytes to provide a relative quantification of each of the two or more different analytes, e.g., relative to each other. In addition, a relative quantitation may be ascertained using a control, or reference, sample as is commonly done in array based assays as well as in quantitative PCR/RT-PCR analyses (described in further detail below).

As noted above, genes/proteins that find use in monitoring a subject for an AR response, i.e., genes/proteins that are differentially expressed or present at different levels in subjects that have had, are currently having or will have an AR episode, are shown ion Tables 1 and 2. Note that for the genes in these tables, detailed information, including precise sequence information, can be determined through the NCBI Entrez Gene database located at the website http(colon)//www(dot)ncbi.nlm.nih(dot)gov. The detailed information for each gene is then obtained by selecting “Gene” and searching for the GeneID No. listed in these tables.

In certain embodiments, additional genes beyond those listed in Tables 1 and 2 may be assayed, and include other genes whose expression level/pattern can be used to determine an AR response as well as genes whose expression level/pattern can be used to evaluate additional transplant characteristics, including but not limited to: a graft tolerant phenotype in a subject, chronic allograft injury (chronic rejection) in blood; immunosuppressive drug toxicity or adverse side effects including drug-induced hypertension; age or body mass index associated genes that correlate with renal pathology or account for differences in recipient age-related graft acceptance; immune tolerance markers in whole blood; genes found in literature surveys with immune modulatory roles that may play a role in transplant outcomes. In addition, other array assay function related genes may be evaluated, e.g., for assessing sample quality (3′- to 5′-bias in probe location), sampling error in biopsy-based studies, cell surface markers, and normalizing genes for calibrating hybridization results (exemplary genes in these categories can be found in U.S. patent application Ser. No. 11/375,681, filed on Mar. 3, 2006, which is incorporated by reference herein in its entirety).

In certain embodiments, additional genes are evaluated to determine whether a subject who has received an allograft has a graft tolerant phenotype, e.g., as described in provisional patent application 61/089,805, filed on Aug. 18, 2008, which is incorporated herein by reference in its entirety. By graft tolerant phenotype is meant that the subject does not reject a graft organ, tissue or cell(s) that has been introduced into/onto the subject. In other words, the subject tolerates or maintains the organ, tissue or cell(s) that has been transplanted to it. A feature of the graft tolerant phenotype detected or identified is that it is a phenotype which occurs without immunosuppressive therapy, i.e., it is present in a subject that is not undergoing immunosuppressive therapy such that immunosuppressive agents are not being administered to the host.

In these embodiments, in addition to evaluating the one or more genes in Tables 1 and 2 to determine an AR response in a subject having an allograft, one or more genes selected from Tables 3 and 4 below are also evaluated to determine whether a subject has a graft tolerant phenotype. Evaluation of any combination of genes in Tables 3 and 4 find use in such embodiments, including evaluating the expression of all genes listed in either or both of Tables 3 and 4.

TABLE 3 List of 28 genes whose expression levels predict 3 distinct subject classes (HD vs. TOL vs. CAN). HD vs. TOL vs. Symbol Entrez Gene ID CAN TP73 7161 Y CIRBP 1153 Y VNN1 8876 Y ATXN1 6310 Y BMP2K 55589 Y KIAA1324 57535 Y PLXNC1 10154 Y CCL3 6348 Y WNK1 65125 Y CCL3L3 414062 Y IL6 3569 Y MYL2 4633 Y KIF15 56992 Y IPO13 9670 Y TMEM117 84216 Y KIAA1609 57707 Y RUNX2 860 Y CBFA2T2 9139 Y TUBB4 10382 Y PDPN 10630 Y PRPF40B 25766 Y NRTN 4902 Y KRT17 3872 Y TSPAN7 7102 Y CLGN 1047 Y CCL2 6347 Y CENPN 55839 Y DOPEY2 9980 Y

TABLE 4 A gene list (24 genes) of significant differentially expressed either up or down among 3-class (HD vs. SP vs. IN) gene Entrez Gene ID Z_JT P value JT HD->SP->IN CD3+ CD4+/CD8+ CD45+/CD45− AGBL3 340351 −3.50419 0.00046 Y ELAVL1 1994 −3.50419 0.00046 Y NUPL1 9818 −3.39854 0.00068 Y KLF5 688 −3.22245 0.00127 Y MYLIP 29116 −3.22245 0.00127 Y ATXN3 4287 −3.15201 0.00162 Y ZDHHC17 23390 −2.97592 0.00292 Y PGM3 5238 −2.87027 0.0041 Y ZMYM2 7750 −2.76462 0.0057 Y FAM110C 642273 −2.58853 0.00964 Y Naïve OR10J3 441911 3.25767 0.00112 Y SSX3 10214 3.25767 0.00112 Y GNG13 51764 3.29289 0.00099 Y KIAA1751 85452 3.3281 0.00087 Y RGMA 56963 3.36332 0.00077 Y IL1F8 27177 3.43376 0.0006 Y GPHA2 170589 3.53941 0.0004 Y MKL2 57496 3.57463 0.00035 Y TRPV1 7442 3.64507 0.00027 Y KRTAP13-1 140258 3.75072 0.00018 Y NOTCH4 4855 3.78594 0.00015 Y OR13C4 138804 3.82116 0.00013 Y KCNE1 3753 3.96203 0.00007 Y

In certain embodiments, the expression profile obtained is a genomic or nucleic acid expression profile, where the amount or level of one or more nucleic acids in the sample is determined, e.g., the nucleic acid transcript of the gene of interest. In these embodiments, the sample that is assayed to generate the expression profile employed in the diagnostic methods is one that is a nucleic acid sample. The nucleic acid sample includes a plurality or population of distinct nucleic acids that includes the expression information of the phenotype determinative genes of interest of the cell or tissue being diagnosed. The nucleic acid may include RNA or DNA nucleic acids, e.g., mRNA, cRNA, cDNA etc., so long as the sample retains the expression information of the host cell or tissue from which it is obtained. The sample may be prepared in a number of different ways, as is known in the art, e.g., by mRNA isolation from a cell, where the isolated mRNA is used as is, amplified, employed to prepare cDNA, cRNA, etc., as is known in the differential expression art. In certain embodiments, the sample is prepared from a cell or tissue harvested from a subject to be diagnosed, e.g., via biopsy of tissue, using standard protocols, where cell types or tissues from which such nucleic acids may be generated include any tissue in which the expression pattern of the to be determined phenotype exists, including, but not limited to, peripheral blood lymphocyte cells, etc., as reviewed above.

The expression profile may be generated from the initial nucleic acid sample using any convenient protocol. While a variety of different manners of generating expression profiles are known, such as those employed in the field of differential gene expression analysis, one representative and convenient type of protocol for generating expression profiles is array-based gene expression profile generation protocols. In certain embodiments, such applications are hybridization assays in which a nucleic acid array that displays “probe” nucleic acids for each of the genes to be assayed/profiled in the profile to be generated is employed. In these assays, a sample of target nucleic acids is first prepared from the initial nucleic acid sample being assayed, where preparation may include labeling of the target nucleic acids with a label, e.g., a member of signal producing system. Following target nucleic acid sample preparation, the sample is contacted with the array under hybridization conditions, whereby complexes are formed between target nucleic acids that are complementary to probe sequences attached to the array surface. The presence of hybridized complexes is then detected, either qualitatively or quantitatively. Specific hybridization technology which may be practiced to generate the expression profiles employed in the subject methods includes the technology described in U.S. Pat. Nos. 5,143,854; 5,288,644; 5,324,633; 5,432,049; 5,470,710; 5,492,806; 5,503,980; 5,510,270; 5,525,464; 5,547,839; 5,580,732; 5,661,028; 5,800,992; the disclosures of which are herein incorporated by reference; as well as WO 95/21265; WO 96/31622; WO 97/10365; WO 97/27317; EP 373 203; and EP 785 280. In these methods, an array of “probe” nucleic acids that includes a probe for each of the phenotype determinative genes whose expression is being assayed is contacted with target nucleic acids as described above. Contact is carried out under hybridization conditions, e.g., stringent hybridization conditions, and unbound nucleic acid is then removed.

The resultant pattern of hybridized nucleic acid provides information regarding expression for each of the genes that have been probed, where the expression information is in terms of whether or not the gene is expressed and, typically, at what level, where the expression data, i.e., expression profile (e.g., in the form of a transcriptosome), may be both qualitative and quantitative.

Alternatively, non-array based methods for quantitating the levels of one or more nucleic acids in a sample may be employed, including quantitative PCR, real-time quantitative PCR, and the like. (For general details concerning real-time PCR see Real-Time PCR: An Essential Guide, K. Edwards et al., eds., Horizon Bioscience, Norwich, U.K. (2004)).

Where the expression profile is a protein expression profile, any convenient protein quantitation protocol may be employed, where the levels of one or more proteins in the assayed sample are determined. Representative methods include, but are not limited to: proteomic arrays, flow cytometry, standard immunoassays (e.g., ELISA assays, western blots), protein activity assays, including multiplex protein activity assays, etc.

Following obtainment of the expression data, or expression profile, from the sample being assayed, the expression profile is analyzed. In certain embodiments, analysis includes comparing the expression profile with a reference or control profile to monitor the subject from which the sample was obtained/derived for an AR episode. The terms “reference” and “control” as used herein mean a standardized pattern of gene expression or levels of expression of certain genes to be used to interpret the expression signature of a given patient with respect to an AR response. The reference or control profile may be a profile that is obtained from a cell/tissue known to have the desired phenotype, e.g., having an AR phenotype, and therefore may be a positive reference or control profile. In addition, the reference/control profile may be from a cell/tissue known to not have the desired phenotype, e.g., an AR negative phenotype (e.g., STA or non-transplanted control), and therefore be a negative reference/control profile.

In certain embodiments, the obtained expression profile is compared to a single reference/control profile to obtain information regarding an AR response. In yet other embodiments, the obtained expression profile is compared to two or more different reference/control profiles to obtain more in depth information regarding the AR phenotype of the assayed cell/tissue. For example, the obtained expression profile may be compared to a positive and negative reference profile to obtain confirmed information regarding whether the cell/tissue has a particular AR phenotype.

The comparison of the obtained expression profile and the one or more reference/control profiles may be performed using any convenient methodology, where a variety of methodologies are known to those of skill in the array art, e.g., by comparing digital images of the expression profiles, by comparing databases of expression data, etc. Patents describing ways of comparing expression profiles include, but are not limited to, U.S. Pat. Nos. 6,308,170 and 6,228,575, the disclosures of which are herein incorporated by reference. Methods of comparing expression profiles are also described above.

The comparison step results in information regarding how similar or dissimilar the obtained expression profile is to the control/reference profile(s), which similarity/dissimilarity information is employed to determine the phenotype of the cell/tissue being assayed and thereby providing a way to monitor AR in the subject. For example, similarity of the obtained gene expression profile with the gene expression profile of a control sample from a subject experiencing an active AR response indicates that the assayed cell/tissue are from a subject experiencing AR. Likewise, similarity of the obtained gene expression profile with the gene expression profile of a control sample from a subject that has not had (or isn't having) an AR episode (e.g., STA) indicates that the assayed cell/tissue are from a subject not experiencing AR.

Depending on the type and nature of the reference/control profile(s) to which the obtained expression profile is compared, the above comparison step yields a variety of different types of information regarding the cell/tissue that is assayed. As such, the above comparison step can yield a positive/negative determination of an ongoing AR response. In many embodiments, the above-obtained information about the cell/tissue being assayed is employed to predict whether a subject will experience an AR response, e.g., within 3 to 6 months from the time of harvesting the sample, and/or to characterize an AR response, e.g., whether an AR response is steroid resistant or sensitive. In certain embodiments, the determination/prediction of AR can be coupled with a determination of additional characteristics of the graft and function thereof. For example, in certain embodiments one can assay for other graft-related pathologies, e.g., chronic rejection (or CAN) and/or drug toxicity (DT) (see, e.g., U.S. patent application Ser. No. 11/375,681, filed on Mar. 3, 2006, which is incorporated by reference herein in its entirety).

The subject methods further find use in pharmacogenomic applications. In these applications, a subject/host/patient is first monitored for AR according to the subject invention, and then treated using a protocol determined, at least in part, on the results of the monitoring. For example, a host may be evaluated for the presence or absence of AR using a protocol such as the diagnostic protocol described in the preceding section. The subject may then be treated using a protocol whose suitability is determined using the results of the monitoring step. For example, where the subject is predicted to have an AR response within the next 3 to 6 months, immunosuppressive therapy can be modulated, e.g., increased or drugs changed, as is known in the art for the treatment/prevention of AR. Likewise, where the subject is predicted to be free of current and near-term AR, the immunosuppressive therapy can be reduced in order to reduce the potential for DT.

In practicing the subject methods, a subject is typically monitored for AR following receipt of a graft or transplant. The subject may be screened once or serially following transplant receipt, e.g., weekly, monthly, bimonthly, half-yearly, yearly, etc. In certain embodiments, the subject is monitored prior to the occurrence of an AR episode. In certain other embodiments, the subject is monitored following the occurrence of an AR episode.

The subject methods may be employed with a variety of different types of transplant subjects. In many embodiments, the subjects are within the class mammalian, including the orders carnivore (e.g., dogs and cats), rodentia (e.g., mice, guinea pigs, and rats), lagomorpha (e.g. rabbits) and primates (e.g., humans, chimpanzees, and monkeys). In certain embodiments, the animals or hosts, i.e., subjects (also referred to herein as patients) are humans.

The methods may be used to monitor a variety of different types of grafts. Grafts of interest include, but are not limited to: transplanted heart, kidney, lung, liver, pancreas, pancreatic islets, brain tissue, stomach, large intestine, small intestine, cornea, skin, trachea, bone, bone marrow, muscle, bladder or parts thereof.

Databases of Expression Profiles of Phenotype Determinative Genes

Also provided are databases of expression profiles of AR responses. Such databases will typically comprise expression profiles of specific tissues from a transplant recipient that are indicative of one or more of: a near-term AR event (within 3 to 6 months), an ongoing AR response, a previous AR response, and a characteristic of an AR response (e.g., steroid resistant/sensitive AR response).

The expression profiles and databases thereof may be provided in a variety of media to facilitate their use (e.g., in a user-accessible/readable format). “Media” refers to a manufacture that contains the expression profile information of the present invention. The databases of the present invention can be recorded on computer readable media, e.g. any medium that can be read and accessed directly by a user employing a computer. Such media include, but are not limited to: magnetic storage media, such as floppy discs, hard disc storage medium, and magnetic tape; optical storage media such as CD-ROM; electrical storage media such as RAM and ROM; and hybrids of these categories such as magnetic/optical storage media. One of skill in the art can readily appreciate how any of the presently known computer readable mediums can be used to create a manufacture comprising a recording of the present database information. “Recorded” refers to a process for storing information on computer readable medium, using any such methods as known in the art. Any convenient data storage structure may be chosen, based on the means used to access the stored information. A variety of data processor programs and formats can be used for storage, e.g. word processing text file, database format, etc. Thus, the subject expression profile databases are accessible by a user, i.e., the database files are saved in a user-readable format (e.g., a computer readable format, where a user controls the computer).

As used herein, “a computer-based system” refers to the hardware means, software means, and data storage means used to analyze the information of the present invention. The minimum hardware of the computer-based systems of the present invention comprises a central processing unit (CPU), input means, output means, and data storage means. A skilled artisan can readily appreciate that any one of the currently available computer-based system are suitable for use in the present invention. The data storage means may comprise any manufacture comprising a recording of the present information as described above, or a memory access means that can access such a manufacture.

A variety of structural formats for the input and output means can be used to input and output the information in the computer-based systems of the present invention, e.g., to and from a user. One format for an output means ranks expression profiles possessing varying degrees of similarity to a reference expression profile. Such presentation provides a skilled artisan (or user) with a ranking of similarities and identifies the degree of similarity contained in the test expression profile to one or more references profile(s).

Reagents, Systems and Kits

Also provided are reagents, systems and kits thereof for practicing one or more of the above-described methods. The subject reagents, systems and kits thereof may vary greatly. Reagents of interest include reagents specifically designed for use in production of the above-described expression profiles of AR phenotype determinative genes, i.e., a gene expression evaluation element made up of one or more reagents. The term system refers to a collection of reagents, however compiled, e.g., by purchasing the collection of reagents from the same or different sources. The term kit refers to a collection of reagents provided, e.g., sold, together.

One type of such reagent is an array of probe nucleic acids in which the phenotype determinative genes of interest are represented. A variety of different array formats are known in the art, with a wide variety of different probe structures, substrate compositions and attachment technologies (e.g., dot blot arrays, microarrays, etc.). Representative array structures of interest include those described in U.S. Pat. Nos. 5,143,854; 5,288,644; 5,324,633; 5,432,049; 5,470,710; 5,492,806; 5,503,980; 5,510,270; 5,525,464; 5,547,839; 5,580,732; 5,661,028; 5,800,992; the disclosures of which are herein incorporated by reference; as well as WO 95/21265; WO 96/31622; WO 97/10365; WO 97/27317; EP 373 203; and EP 785 280.

In certain embodiments, the arrays include probes for at least 1 of the genes listed in Table 1. As such, probes for any combination of genes in Table 1 may be employed. Therefore, in certain embodiments, the number of genes that are from Table 1 that are represented on the array is at least 2, at least 3, at least 4, at least 5, at least 8 or more, including all of the genes listed in Table 1. The subject arrays may include only those genes that are listed in Table 1 or they may include additional genes that are not listed in Table 1, such as probes for genes whose expression pattern can be used to evaluate additional transplant characteristics, including but not limited to: chronic allograft injury (chronic rejection) in blood; immunosuppressive drug toxicity or adverse side effects including drug-induced hypertension; age or body mass index associated genes that correlate with renal pathology or account for differences in recipient age-related graft acceptance; immune tolerance markers in whole blood; genes found in literature surveys with immune modulatory roles that may play a role in transplant outcomes; as well as other array assay function related genes, e.g., for assessing sample quality (3′- to 5′-bias in probe location), sampling error in biopsy-based studies, cell surface markers, and normalizing genes for calibrating hybridization results; and the like. Where the subject arrays include probes for such additional genes, in certain embodiments the number % of additional genes that are represented and are not directly or indirectly related to transplantation does not exceed about 50%, usually does not exceed about 25%. In certain embodiments where additional genes are included, a great majority of genes in the collection are transplant characterization genes, where by great majority is meant at least about 75%, usually at least about 80% and sometimes at least about 85, 90, 95% or higher, including embodiments where 100% of the genes in the collection are phenotype determinative genes. Transplant characterization genes are genes whose expression can be employed to characterize transplant function in some manner, e.g., presence of rejection, etc.

Another type of reagent that is specifically tailored for generating expression profiles of phenotype determinative genes is a collection of gene specific primers that is designed to selectively amplify such genes (e.g., using a PCR-based technique, e.g., real-time RT-PCR). Gene specific primers and methods for using the same are described in U.S. Pat. No. 5,994,076, the disclosure of which is herein incorporated by reference. Of particular interest are collections of gene specific primers that have primers for at least 1 of the genes listed in one Table 1, often a plurality of these genes, e.g., at least 2, 4, 8 or more. In certain embodiments, all of genes that are from Table 1 have primers in the collection. The subject gene specific primer collections may include only those genes that are listed in Table 1, or they may include primers for additional genes that are not listed in Table 1, such as probes for genes whose expression pattern can be used to evaluate additional transplant characteristics, including but not limited to: chronic allograft injury (chronic rejection) in blood; immunosuppressive drug toxicity or adverse side effects including drug-induced hypertension; age or body mass index associated genes that correlate with renal pathology or account for differences in recipient age-related graft acceptance; immune tolerance markers in whole blood; genes found in literature surveys with immune modulatory roles that may play a role in transplant outcomes; as well as other array assay function related genes, e.g., for assessing sample quality (3′- to 5′-bias in probe location), sampling error in biopsy-based studies, cell surface markers, and normalizing genes for calibrating hybridization results; and the like. Where the subject arrays include probes for such additional genes, in certain embodiments the number % of additional genes that are represented and are not directly or indirectly related to transplantation does not exceed about 50%, usually does not exceed about 25%. In certain embodiments where additional genes are included, a great majority of genes in the collection are transplant characterization genes, where by great majority is meant at least about 75%, usually at least about 80% and sometimes at least about 85, 90, 95% or higher, including embodiments where 100% of the genes in the collection are phenotype determinative genes.

The systems and kits of the subject invention may include the above-described arrays and/or gene specific primer collections. The systems and kits may further include one or more additional reagents employed in the various methods, such as primers for generating target nucleic acids, dNTPs and/or rNTPs, which may be either premixed or separate, one or more uniquely labeled dNTPs and/or rNTPs, such as biotinylated or Cy3 or Cy5 tagged dNTPs, gold or silver particles with different scattering spectra, or other post synthesis labeling reagent, such as chemically active derivatives of fluorescent dyes, enzymes, such as reverse transcriptases, DNA polymerases, RNA polymerases, and the like, various buffer mediums, e.g. hybridization and washing buffers, prefabricated probe arrays, labeled probe purification reagents and components, like spin columns, etc., signal generation and detection reagents, e.g. streptavidin-alkaline phosphatase conjugate, chemifluorescent or chemiluminescent substrate, and the like.

The subject systems and kits may also include a phenotype determination element, which element is, in many embodiments, a reference or control expression profile that can be employed, e.g., by a suitable computing means, to make an AR phenotype determination based on an “input” expression profile, e.g., that has been determined with the above described gene expression evaluation element. Representative phenotype determination elements include databases of expression profiles, e.g., reference or control profiles, as described above.

In addition to the above components, the subject kits will further include instructions for practicing the subject methods. These instructions may be present in the subject kits in a variety of forms, one or more of which may be present in the kit. One form in which these instructions may be present is as printed information on a suitable medium or substrate, e.g., a piece or pieces of paper on which the information is printed, in the packaging of the kit, in a package insert, etc. Yet another means would be a computer readable medium, e.g., diskette, CD, etc., on which the information has been recorded. Yet another means that may be present is a website address which may be used via the internet to access the information at a removed site. Any convenient means may be present in the kits.

The following examples are offered by way of illustration and not by way of limitation.

EXPERIMENTAL Example I Methods

Study Design

Patient demographics and clinical data were collected for all recipients providing peripheral blood samples in this study. The mean±standard deviation data for the different study sets (Verification Set, Validation Set 1 and Validation Set 2) is provided in FIG. 11 and Table 1 (below). Within each sample set, the AR and STA phenotypes have matched demographics. To interrogate the performance of blood-based genes for AR diagnosis, irrespective of time post-transplantation, sample study sets were deliberately chosen to have different mean times post-transplantation.

TABLE 1 qPCR Microarray_Affy Verification Set Clinical AR STA p AR STA p Characteristics (n = 31) (n = 28) value (n = 17) (n = 17) value Recipients Gender, % females 34.3% 51.4% 0.1 34.6% 50% 0.3 Mean age, year 12.0 ± 5.1  10.7 ± 6.2  0.2 11.8 ± 4.7  7.8 ± 6.8 0.1 Immunosuppression, % SF 45.7% 54.1% 0.5 46.2% 50.0% 0.8 HLA match 1.2 ± 1.7 1.4 ± 1.1 0.7 0.8 ± 0.8 0.7 ± 1.1 0.9 Post-txp time, month 34.8 ± 40.7 23.4 ± 31.2 0.1 31.1 ± 39.4   14 ± 22.4 0.1 Donor Donor Source % LRD 60.0% 62.2% 0.9 57.7% 85.7% 0.1 Gender, % females 53.1% 48.7% 0.7 50.0% 42.9% 0.7 Mean age, year 34.9 ± 11.3 30.8 ± 11.4 0.1 33.3 ± 10.9 30.6 ± 10.4 0.5 qPCR Validation Set 1 Validation Set 2 Clinical AR STA p AR STA p Characteristics (n = 21) (n = 21) value (n = 32) (n = 32) value Recipients Gender, % females 31.6% 42.9% 0.5 33.3% 43.8% 0.5 Mean age, year 11.2 ± 5.9  13.5 ± 5   0.2 12.5 ± 6.2  13.9 ± 6.3  0.5 Immunosuppression, % SF 63.1% 57.1% 0.7 66.7% 68.8% 0.9 HLA match 1.6 ± 1.6 0.8 ± 1.0 0.1 1.6 ± 1.3 1.8 ± 0.9 0.2 Post-txp time, month 12.4 ± 17.9 10.6 ± 9   0.7 7.9 ± 5.2 6.5 ± 4.4 0.4 Donor Donor Source % LRD 31.6% 61.9% 0.1 55.6% 75.0% 0.2 Gender, % females 31.6% 57.1% 0.1 44.4% 50.0% 0.7 Mean age, year 26.4 ± 11.8 31.1 ± 10.9 0.2 22.1 ± 11.7   37 ± 10.8 0.2 Values are means ± SD (Standard Deviation). AR, Acute rejection; STA, stable; SF, steriod free; txp, transplant; LRD, living related donor

This study was undertaken as a sum of 3 projects (FIG. 1): In the first project, we performed cross-sectional microarray analysis of 122 peripheral blood samples collected either as 2.5 ml of whole blood samples (n=106) or as 5 ml of blood (n=16) for isolation of peripheral blood leukocytes (PBL)⁷. Blood samples were collected at a single time point, timed with a biopsy where acute rejection was either confirmed¹⁰, as present (60 AR samples) or absent (62 STA samples). Sample preparation for array hybridization was done following the manufacturers' recommended protocols.

The samples were hybridized to one of 3 microarray platforms: Affymetrix (n=75, the 54K HGU133Plus2 Array), Agilent (n=26, 44K Whole Human Genome oligo Array) or the Lymphochip cDNA array (n=21, 32K cDNA clones). Appropriate probe detection and normalization methods were applied for each microarray platform. Details of sample preparation, array hybridization and analyses¹¹ ¹² ¹³ ¹⁴ are provided in Supplemental Methods (see below). All microarray data files are available in GEO under Accession No. GSE14067.

In the second project, Q-PCR analysis was performed for verification and validation on selected genes from the array studies in project 1. Expression of messenger RNA for the selected candidate genes was measured in each blood sample by Q-PCR using TaqMan Gene Expression Assay products (Applied Biosystems, Foster City, Calif.). 32 genes were selected for Q-PCR verification (see Supplementary Methods/Candidate Gene Selection Criteria for Q-PCR) in 34 samples (17 AR, 17 STA; Verification Set) previously used for microarray, left with adequate quantity and quality residual RNA. Of the 32 genes, the 10 most significant genes were selected for further validation on 42 independent samples (21 AR, 21 STA; Validation Set 1) and an additional 64 independent samples (32 AR, 32 STA; Validation Set 2). We applied various classification models¹⁵ including linear discriminant analysis, 5 different Bayesian models, multilayer perceptron and logistic regression (multinomial and linear) models to the 10 gene-set expression values obtained from Validation Set 1. Each of these models was built using Validation Set 1 and tested on the Validation Set 2 for evaluating their sensitivities, specificities and accuracies. A multinomial logistic regression model using 5 out of the 10 gene-set yielded the best overall performance. The primer list for each of the 10 genes, their corresponding ABI probe IDs, the expression of the 10 most significant genes by Q-PCR, and the comparisons of all classification models used, are all provided in FIGS. 7 to 10.

In the third project, we applied the 5-gene AR diagnosis model for AR prediction modeling on 40 patients who had serially collected peripheral blood samples, collected either before (preAR; n=27 samples) and/or after (postAR; n=30 samples) the AR episode. Longitudinal 5-gene Q-PCR analysis was done on paired preAR blood samples, collected 1 to 6 months prior to the AR episode, to estimate the prediction of AR injury, even at the time of unperturbed graft function. These were the only samples in the study that were not paired with graft biopsies. A similar analysis was done on postAR samples collected 1 to 6 months after treatment of the AR episode (paired with graft biopsies performed to follow-up for AR resolution) to interrogate the effect of intensified immunosuppression on the 5 gene-set signature (FIG. 1).

A first logistic regression model using the 5-genes set is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {10.0395 + {({1.8364*{DUSP}\; 1})} + {({{- 1.9182}*{MAPK}\; 9})} + {({{- 1.8939}*{NKTR}})} +} \\ {{({4.9901*{PBEF}\; 1})} + {({2.124*{PSEN}\; 1})}} \end{matrix}\;}}{1 + ^{\begin{matrix} {10.0395 + {({1.8364*{DUSP}\; 1})} + {({{- 1.9182}*{MAPK}\; 9})} +} \\ {{({{- 1.8939}*{NKTR}})} + {({4.9901*{PBEF}\; 1})} + {({2.124*{PSEN}\; 1})}} \end{matrix}}}$

A second logistic regression model using a minimal 5-gene set is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {{- 18.6137} + {({3.8327*{DUSP}\; 1})} + {({{- 5.1524}*{MAPK}\; 9})} +} \\ {{({{- 3.6056}*{NKTR}})} + {({8.3652*{PBEF}\; 1})} + {({3.353*{PSEN}\; 1})}} \end{matrix}}}{1 + ^{\begin{matrix} {{- 18.6137} + {({3.8327*{DUSP}\; 1})} + {({{- 5.1524}*{MAPK}\; 9})} +} \\ {{({{- 3.6056}*{NKTR}})} + {({8.3652*{PBEF}\; 1})} + {({3.353*{PSEN}\; 1})}} \end{matrix}}}$

The logistic regression model using a minimal 5-gene set based on the Stanford samples is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {{- 17.2226} + {({1.0735*{IFNGR}\; 1})} + {({{- 7.7795}*{MAPK}\; 9})} + {({4.2314*{PBEF}\; 1})} +} \\ {{({3.7546*{PSEN}\; 1})} + {({2.9103*{RNF}\; 130})}} \end{matrix}}}{1 + ^{\begin{matrix} {{- 17.2226} + {({1.0735*{IFNGR}\; 1})} + {({{- 7.7795}*{MAPK}\; 9})} + {({4.2314*{PBEF}\; 1})} +} \\ {{({3.7546*{PSEN}\; 1})} + {({2.9103*{RNF}\; 130})}} \end{matrix}}}$

In each equation above, θ is the predicted probability of a given sample belonging to one of the classes (AR or STA). The threshold θ=0.5 was selected for the best sensitivity and specificity, based on the Receiver Operating Characteristic (ROC) curve, to determine whether the predictive class was an AR or STA phenotype. The model was tested in Validation Set 1 and then re-validated in Validation Set 2, for true AR classification. The accuracy of the model was assessed by evaluating the sensitivity, specificity, positive prediction value (PPV) and negative prediction value (NPV) on the test set (Validation set 2). Finally, after confirming the robustness of the model for AR diagnosis, the model was used to predict the probability of AR on paired blood samples, sequentially collected before and after AR.

Statistical and Bioinformatics Analyses

Student t-test was used for establishing the significance for gene expression differences between AR and STA samples by Q-PCR. All differences were considered statistically significant for false discovery rate (FDR)<10%. Data was collected on the following 18 demographic and clinical variables on all samples processed for Q-PCR validation studies: post-transplant time, recipient age, recipient gender, donor gender, donor source, donor age, cytomegalovirus/CMV, Epstein-Barr virus/EBV and BK virus/BKV viremia at sampling, concomitant bacterial infection, presence of donor specific antibodies (DSA) or panel reactive antibodies (PRA), usage of induction therapy, maintenance steroid usage, use of different calcineurin inhibitor drugs (FK506 or cyclosporine A/CSA), use of different anti-metabolites (mycophenolate mofetil/MMF or azathioprine/AZA), the white blood cell (WBC) count and the hematocrit (HCT) levels at the time of blood sampling (FIG. 11). Data was provided for SNSO1 samples by PPD (David Ikle). Pearson correlation and t-test were used, as appropriate, to determine any influence of the 18 parameters as potential confounders for expression of the 5-genes set (FIG. 12). As most samples had paired biopsies for evaluation, samples were divided into 2 groups according to Banff score grades (1=Banff scores≦1A; 2=Banff score≧1B); C4d⁺ was determined if ≧0% of peritubular capillaries (PTC) had C4d deposits (FIG. 14).

Biological relevance of the cross-platform 525 gene set was evaluated using the PANTHER classification system (http(colon)//www(dot)pantherdb(dot)org) and Ingenuity Pathway Analysis (IPA; http(colon)//www(dot)ingenuity(dot)com) (see FIGS. 16 and 17). More information is provided in Supplemental Methods.

Results

Significant Genes for AR in Peripheral Blood by Cross-Platform Microarray Analysis

Microarray data generated from 3 array platforms were cross-compared, by mapping common and overlapping transcripts to Human Gene Organization (HUGO) gene names. Because each array platform uses different sets of genes, which are represented using different probe IDs, we used AILUN (http(colon)//ailun(dot)stanford(dot)edu)¹² to re-annotate the probe identifiers (IDs) with the current HUGO gene names, which identified 10,357 common genes that were shared across all 3 platforms. We applied supervised, two-class unpaired Significance Analysis of Microarray (SAM¹¹) analyses with a common significance threshold of false discovery rate <10% to each of the 3 re-annotated array lists. We identified 9710 genes on the Affymetrix arrays, 8642 genes on the Agilent arrays and 2805 genes on the cDNA arrays as significantly differentially expressed in AR. Despite the large number of significant genes on each platform, only 525 genes were differentially expressed on all 3 platforms (FIG. 5), irrespective of the blood sample collection procedure.

Q-PCR Verification and Validation (Validation Set 1)

32 genes were selected for initial Q-PCR verification in 34 samples (17 AR, 17 STA) that had also been used on the arrays studies. From this data, the 10 most significant genes (q<10%) were further validated in 42 independent blood samples (Validation Set 1; 21 AR, 21 STA). The expression fold changes of these 10 most significant genes were in agreement with array data, with same direction fold changes and higher signal by Q-PCR than by microarray (FIG. 1), and separation of the samples into AR or STA phenotypes by unsupervised clustering of expression data for the 10-gene set by Q-PCR (FIG. 2). Individual expression data for each gene and each sample is shown in FIG. 6B. For the 3 genes chosen from the AR literature, Q-PCR expression of FOXP3¹⁶ and GRZYB¹⁷ ¹⁸ were increased in AR compared to STA blood samples (p=0.001, p=0.03 respectively) (FIG. 7), even though these genes were not among the significant gene list from Affymetrix microarrays, and both PFN1¹⁹ and GRZYB were conversely down regulated in AR on the cDNA and Agilent platform respectively, with no statistical difference by Q-PCR between AR and STA samples. Hence none of these 3 genes were used for further cross-validation studies.

Cross-Validation (Validation Set 2) and Prediction Modeling for AR Classification

A multinomial logistic regression model, which demonstrated the best performance of the 10 classification models tested (details in FIG. 10), resulted in a minimal 5 gene model for AR prediction with high confidence shown by ROC score=95.2%, when tested on Q-PCR data from the 64 independent samples in Validation Set 2. The 5 gene-set expression shows excellent AR and STA class separation (FIG. 3), with a high AR prediction score with 100% sensitivity, 93.6% specificity, 94.1% positive predictive value (PPV) and 100% negative predictive value (NPV). The only misclassification related to 2 STA samples. Both misclassified STA samples developed clinical AR at 4 and 7 months after the initial sampling for this study, suggesting that the AR prediction of the 5 gene-test maybe more robust than estimated by the prediction model.

Evaluating for Demographic, Clinical and Biopsy Confounders for the 5 Gene-Set for AR Diagnosis

To examine if any baseline, or at sampling, demographic, clinical or immunosuppression confounders may drive the segregation of the 5 gene-set prediction score for AR, Pearson correlation coefficients calculated for data collected on 18 different confounders (FIG. 11), showed no confounding effects on the expression of the 5 gene-set (maximum |r|=0.27) (FIGS. 12 and 13). By t-test, DUSP1 and PSEN1 had higher expression (p<0.01) in patients with positive DSA. Similarly, we also examined if the grade or type of rejection had an influence on expression of the 5 gene-set (FIG. 14). No correlation was found between AR probability prediction by the 5 gene-set and Banff scores (r=−0.13), humoral versus cellular rejection (r=0.07), and presence or absence of peri-tubular C4d+ (r=0.2) (FIG. 14). These data support that the coordinated expression the 5 gene-set in peripheral blood can diagnose AR with high confidence, irrespective of the differences in patient characteristics and commonly used immunosuppression choices, rejection timing or the underlying mechanism of rejection injury.

Prediction of AR Injury Prior to Graft Dysfunction and Clinically Indicated Biopsies

We applied the AR regression model built on cross-sectional samples, to 27 preAR blood samples (sequentially collected 1 to 6 month prior to a biopsy proven AR episode, prior to any clinical graft dysfunction), as part of serial blood sampling protocol post-transplantation. Additionally, 70 serial blood samples from 58 unique STA patients, who did not have any AR, collected at the time of protocol biopsies at 1, 3, 6, 12 and 24 months post-transplantation were also profiled by Q-PCR for the 5 gene-set for AR prediction. As shown in FIG. 4, the 5-gene signature from STA samples that were collected at various time points post-transplantation remains “quiescent”, with low AR prediction scores and with little change in expression scores over time post-transplantation (FIG. 4A). However, in the AR patient group (FIG. 4B; preAR, AR and postAR patients), the expression level for the 5 gene-set is at a significantly (p<0.001) higher level in rejector patients (>50% AR prediction scores in all AR patients, versus <20% in all STA patients), even 6 months prior to an AR episode. This finding suggests that patients who undergo AR episodes after transplantation have an immune activation profile in their peripheral blood that can be measured quantitatively well before a clinically relevant and biopsy documented AR episode. It is noteworthy that the AR prediction score for blood samples examined within 3 months of the rejection episode is equivalent to the sample collected at the time of AR (FIG. 4B), providing a potential window for immunosuppression intensification to mitigate or treat early AR and/or instigate an earlier biopsy. When immunosuppression treatment is intensified at AR (steroid bolus or thymoglobulin in this study), there is a rapid decrement of the 5-gene signature over the next 3 months regardless of AR treatment drug (FIG. 4B). FIG. 15 provides mean AR probabilities along with standard errors and calculated creatinine clearances²⁰ for each group of samples.

Discussion

The transplant biopsy is the current gold standard for AR diagnosis but is limited by its invasive nature, difficulty in obtaining multiple samples, and often identifying AR only in the presence of advanced tissue injury and clinical graft dysfunction. There is no reliable blood test currently available to diagnose and predict AR, as many clinical confounders, including over-representation of globin genes in whole blood samples, drive expression changes⁷, and in the past, robust blood-based AR biomarker identification has been hampered by small sample size, experimental variance⁷, concomitant infections and the underlying molecular heterogeneity of tissue rejection injury¹.

As demonstrated above, we have identified a highly stringent gene-set whose expression evaluation in a transplant patient sample (e.g., peripheral blood) can be used to monitor an AR response. This gene set evaluation method proves accurate across multiple different confounders, such as differences in sample preparation, recipient and donor demographics, recipient age, post-transplant time, immunosuppression usage, or concurrent infection. We have validated a minimal set of 5 genes for non-invasive AR diagnosis in 152 independent blood samples. The expression of the 5-gene set is sufficient for prediction of biopsy-proven AR with high confidence levels even 3 months prior to clinical graft dysfunction. Additionally, the expression of the 5-gene set is differentially regulated in rejectors, compared to non-rejectors, up to 7 months prior to a biopsy proven AR episode, indicating that differential regulation of the 5 gene-set in blood can be used as a yardstick for immune activation and risk assessment for future AR. As the 5 gene signature “normalizes” to levels closer to those seen in non-rejector STA patients with intensification of immunosuppression, serial measurement of this gene-set can serve as an immune “barometer” to define an individual patient's immune activation threshold for prediction of AR many months before clinical graft dysfunction. This is an important and critical step towards individualizing immunosuppressive regimens in a transplant recipient.

REFERENCES

1. Sarwal M, Chua M S, Kambham N, et al. Molecular heterogeneity in acute renal allograft rejection identified by DNA microarray profiling. N Engl J Med 2003; 349(2):125-38.

2. Flechner S M, Kurian S M, Solez K, et al. De novo kidney transplantation without use of calcineurin inhibitors preserves renal structure and function at two years. Am J Transplant 2004; 4(11):1776-85.

3. Kurian S M, Flechner S M, Kaouk J, et al. Laparoscopic donor nephrectomy gene expression profiling reveals upregulation of stress and ischemia associated genes compared to control kidneys. Transplantation 2005; 80(8):1067-71.

4. Morgun A, Shulzhenko N, Perez-Diez A, et al. Molecular profiling improves diagnoses of rejection and infection in transplanted organs. Circ Res 2006; 98(12):e74-83.

5. Hotchkiss H, Chu T T, Hancock W W, et al. Differential expression of profibrotic and growth factors in chronic allograft nephropathy. Transplantation 2006; 81(3):342-9.

6. Park W, Griffin M, Grande J P, Cosio F, Stegall M D. Molecular evidence of injury and inflammation in normal and fibrotic renal allografts one year posttransplant. Transplantation 2007; 83(11):1466-76.

7. Li L, Ying L, Naesens M, et al. Interference of globin genes with biomarker discovery for allograft rejection in peripheral blood samples. Physiol Genomics 2008; 32(2):190-7.

8. Deng M C, Eisen H J, Mehra M R, et al. Noninvasive discrimination of rejection in cardiac allograft recipients using gene expression profiling. Am J Transplant 2006; 6(1):150-60.

9. Gwinner W. Renal transplant rejection markers. World journal of urology 2007; 25(5):445-55.

10. Racusen L C, Solez K, Colvin R B, et al. The Banff 97 working classification of renal allograft pathology. Kidney Int 1999; 55(2):713-23.

11. Efron B, Tibshirani R. Empirical bayes methods and false discovery rates for microarrays. Genet Epidemiol 2002; 23(1):70-86.

12. Chen R, Li L, Butte A. AILUN: reannotating gene expression data automatically. Nature Methods 2007; 4(11):In Press.

13. Perez-Iratxeta C, Andrade M A. Inconsistencies over time in 5% of NetAffx probe-to-gene annotations. BMC Bioinformatics 2005; 6:183.

14. Harbig J, Sprinkle R, Enkemann S A. A sequence-based identification of the genes detected by probesets on the Affymetrix U133 plus 2.0 array. Nucleic Acids Res 2005; 33(3):e31.

15. Witten I H, Frank E. Data Mining: Practical machine learning tools and techniques. 2nd ed. San Francisco: Morgan Kaufmann; 2005.

16. Li B, Hartono C, Ding R, et al. Noninvasive diagnosis of renal-allograft rejection by measurement of messenger RNA for perforin and granzyme B in urine. N Engl J Med 2001; 344(13):947-54.

17. Dugre F J, Gaudreau S, Belles-Isles M, Houde I, Roy R. Cytokine and cytotoxic molecule gene expression determined in peripheral blood mononuclear cells in the diagnosis of acute renal rejection. Transplantation 2000; 70(7):1074-80.

18. Simon T, Opelz G, Wiesel M, Ott R C, Susal C. Serial peripheral blood perforin and granzyme B gene expression measurements for prediction of acute rejection in kidney graft recipients. Am J Transplant 2003; 3(9):1121-7.

19. Vasconcellos L M, Schachter A D, Zheng X X, et al. Cytotoxic lymphocyte gene expression in peripheral blood leukocytes correlates with rejecting renal allografts. Transplantation 1998; 66(5):562-6.

20. Schwartz G J, Haycock G B, Edelmann C M, Jr., Spitzer A. A simple estimate of glomerular filtration rate in children derived from body length and plasma creatinine. Pediatrics 1976; 58(2):259-63.

Supplemental Methods

Patients and Samples

274 unique peripheral blood samples from 274 pediatric and young adult kidney allograft recipients were collected at Stanford University (n=193) and other SNS01 centers (n=81) (from 2002 to 2007). All blood samples were obtained at the time of matched clinically indicated or protocol graft biopsies which were used to categorize the sample as acute rejection (AR; n=109)¹, or stable, if there was absence of AR or other substantive pathology (STA, n=108). The diagnosis of acute rejection was confirmed on re-evaluation of all samples by a single pathologist (NK), before their inclusion in this study. AR samples were obtained prior to treatment intensification for rejection. An additional 57 sequential samples were also collected before AR (preAR, n=27) and after AR (postAR, n=30), and serial post-transplant blood samples were also collected from a subset of STA patients (n=70). Written informed consent was obtained from all the subjects. The study was approved by the Stanford Institutional Review Board and by each SNSO1 center Institutional Review Board.

Sample Collection, RNA Extraction, Microarray Target Preparation and Hybridization

Peripheral blood samples were collected in the form of either whole blood (n=106) or peripheral blood leukocytes (PBL; n=16). All SNSO1 samples were collected as whole blood; the samples from Stanford were collected as PBL between 2002-2004, and as either whole blood or PBL, between 2004-2007.

For total RNA extracted from whole blood, 2.5 ml of peripheral blood was collected in PAXgene™ Blood RNA Tubes (PreAnalytiX, Qiagen) and processed using the PAXgene Blood RNA Kit (PreAnalytiX, Qiagen). For RNA extracted from leukocytes, 5-10 ml whole blood was collected in sodium heparin tubes. Leukocytes were recovered by initial separation of the white blood cells by centrifugation and lysis of erythrocytes². RNA was extracted using RNeasy Midi Kit® (Qiagen Inc, Valencia, Calif.). RNA concentration was measured using NanoDrop® ND-1000 (NanoDrop Technologies, Wilmington, Del.) and the integrity of RNA was assessed with the Agilent 2100 Bioanalyzer using RNA Nano Chips (Agilent Technologies, Santa Clara, Calif.). RNA was stored in −80° C. until further use for microarray and Q-PCR.

Sample preparations for Affymetrix and Agilent array hybridizations were done using the T7 protocol (T7), which employs one-cycle cRNA synthesis by using T7-Oligo (dT) Primer. For Affymetrix arrays, a total of 2 μg total RNA, was used for the synthesis of cRNA as the final product following manufactures instructions. cRNA were hybridized on GeneChip® Human Genome U133 Plus 2.0 Arrays (Affymetrix Inc. Santa Clara, Calif.). 100 ng of total RNA was used for Agilent array hybridization following the manufacturer's instruction. Amplified and purified cRNA were incorporated with of Cy5 dye (test sample) or Cy3-dye (common reference as same as used in cDNA array), the mixture of two labeled samples was applied onto either Agilent Whole Human Genome Oligo Microarray 44 k or 4×44 k. cDNA microarrays³, containing 32,000 cDNA clones (12,400 known unique genes), were processed using established protocols, using 2 ug RNA in each channel against a common reference (RNA pool from multiple tumor lines). Hybridized Affymetrix and cDNA microarrays were scanned using GenePix 4000 (Axon Instruments, Union City, Calif.) and fluorescent images were analyzed with the GenePix Pro software package. Agilent slides were scanned by using an Agilent DNA microarray scanner.

Gene Expression Profile Analyses

Different criteria were used for probe detection for different array platforms. For the Affymetrix HG133plus2.0 array, original array data were Quartile-Quartile normalized using dChip software (DNA-Chip Analyzer, www(dot)dchip(dot)org) after probe level intensities and quality were calculated and measured including median intensity, percent of probe set outliers, and percent of single probe outliers. For the Agilent 44 k, a cut off for absolute value of log₂red channel/green channel >0.5 for at least one array was applied; data was normalized in GeneSpring. For the Lymphochip cDNA array, a mean of channel 1 intensity/media background intensity >1.5, and/or normalized mean of channel 2 intensity/media background intensity) >1.5 were used. Data was normalized in the Stanford Microarray Database (Stanford). All Gene expression values were transformed to log2 for further analysis. Unique genes across all platforms were mapped based on mapping of transcripts (probe sets) to Human Gene Organization (HUGO) gene names.

Quantitative Real-Time PCR (Q-PCR)

The probes used for Q-PCR were selected from the Refseq ID for each gene from Applied Biosystems (ABI). The list of probes used is shown in FIG. 6A. cDNA was first transcribed from 500 ng total RNA in 11 ul, and denatured in 1 ul 10 mM dNTPs and 1 ul 300 ng/ul random hexamers at 65° C. for 15 min. 7 ul of the Q-PCR mixture [4 ul 5×binding buffer, 1 ul 0.1 M DTT, 1 ul SuperScript III RT(Invitrogen), 1 ul RNase OUT (Invitrogen)] was added to the denatured RNA tube in a final 20 ul. RT reaction was carried out in a PCR thermal cycler at 50° C. for 60 min, 72° C. for 15 min following by addition of 1 ul RNase H for 30 min at 37° C. to digest away remaining RNA template. Q-PCR was carried out following manufacture's protocol with modification. cDNA was diluted to 1.25 ng/ul. Quantitative RT-PCR reactions were carried out in 384-well plates in a 10 μl volume containing 4 μl of cDNA (5 ng), 1 μl of TaqMAN probe and 5 μl 2× gene expression master mix buffer (ABI, Foster City, Calif.). The Q-PCR reactions were run in an ABI 7700 sequence detector (Applied Biosystems, Foster City, Calif.) under cycle conditions (10 min at 95° C., 40 cycles of 15 s 95° C., 30 s at 60° C.) recommended by manufacturer. The relative amount of RNA expression was calculated using a comparative C_(T) method. Expression values were normalized to 18S ribosomal RNA and control samples (Stratagene Universal RNA) on the plate.

Study Design

This study was undertaken as a sum of 3 projects (FIG. 1): In the first project, we performed cross-sectional microarray analysis of 122 peripheral blood samples collected either as 2.5 ml of whole blood samples (n=106) in PAXgene™ Blood RNA Tubes (PreAnalytiX, Qiagen), or as 5 ml of blood in sodium heparin tubes (n=16) for isolation of peripheral blood leukocytes (PBL). For this analysis, blood samples were collected at a single time point, timed with a biopsy where acute rejection was either confirmed as present (60 AR samples) or absent (62 STA samples). The samples were hybridized to one of 3 microarray platforms: Affymetrix (n=75, the 54K HGU133Plus2 Array), Agilent (n=26, 44K Whole Human Genome oligo Array), or the Lymphochip cDNA array (n=21, 32K cDNA clones). Microarray data generated from 3 array platforms were cross-compared, by mapping common and overlapping transcripts to Human Gene Organization (HUGO) gene names. Significant gene lists on each platform were identified using Significant Analysis of Microarray (SAM,⁴) with a common significance threshold of a false discovery rate (FDR) of <10%. We identified 525 common genes that were significantly differentially expressed in AR on all 3 platforms.

In the second project, 32 genes were selected for Q-PCR verification (see section Candidate Gene Selection Criteria) in 34 samples previously used in microarray (Verification Set). Of the 32 genes tested by Q-PCR, 10 most significant genes were selected for further validation on 42 independent samples (Validation Set 1) and an additional 64 independent samples (32 AR, 32 STA; Validation Set 2). We applied various classification models including linear discriminant analysis, 5 different Bayesian models, multilayer perceptron and logistic regression (multinomial and linear) models to the 10 gene-set in the Validation Set 1. Each of these models was applied to an independent set of 64 samples (Validation Set 2) for evaluating their sensitivities, specificities and accuracies. A multinomial logistic regression model using 5 out of the 10 gene-set yielded the best overall performance.

In the third project, we applied the 5-gene AR diagnosis model, to predict AR in 30 patients who had serially collected peripheral blood samples, collected before (n=27 samples) and after (n=30 samples) the AR episode. Longitudinal 5-gene Q-PCR analysis was done on paired blood samples collected 1 to 6 months prior to AR episode to estimate the prediction of AR injury, even at the time of unperturbed graft function. These were the only samples in the study that were not paired with graft biopsies. A similar analysis was done on samples collected 1 to 6 months after treatment of the AR episode (paired with graft biopsies performed to follow-up for AR resolution) to interrogate the effect of intensified immunosuppression on the 5 gene-set signature (FIG. 1).

Project 1—Discovery of AR Specific Genes Across Multiple Microarray Platforms

Gene Expression Profiling of Peripheral Blood Samples Across 3 Microarray Platforms

Out of the 274 unique peripheral blood samples, 103 unique samples (47 AR, 56 STA) were hybridized on 122 microarrays (60 AR, 62 STA) from 3 different platforms. Out of the 103 peripheral blood samples, 75 were hybridized on Affymetrix HG U133 plus 2.0 microarrays, 26 on Agilent microarrays and 21 on cDNA microarrays. 17 samples were hybridized on 2 or more platforms for quality control purposes. For the samples that were hybridized on more than 1 platform, representative cross-platform correlations were as follows, r=0.6 between cDNA and Agilent, r=0.8 between Affymetrix and Agilent (see cross-reference of samples and assays shown in FIG. 18).

Significance Analysis of Microarrays

We applied significance analysis of microarray (SAM⁴) to identify differentially expressed genes for AR on all 3 platforms. All gene expression values were transformed to log₂ before SAM analysis. False discovery rate (FDR)<10% was used as a threshold on all platforms. Using this threshold we identified 5285 genes on the Affymetrix array, 8642 genes on the Agilent array, and 2805 genes on the cDNA microarrays as differentially expressed genes of significance in AR (FIG. 5). Because each array platform uses different sets of genes, which are represented using different probe IDs, we used AILUN (http://ailun.stanford.edu)⁵ to re-annotate the probe identifiers (IDs) with the current Entrez Gene IDs. Using the Entrez Gene IDs, we identified 525 genes that were differentially expressed on all 3 platforms (FIG. 5). At the same FDR, the cDNA platform seems to have the least number of significant AR genes, when compared to the Affymetrix and Agilent platforms (data not shown).

While these studies show that cross platform array comparisons of a carefully selected phenotypes, is a powerful approach to identify gene specific signatures, a portion of the overlapping genes across all platforms show a converse fold change across the 3 platforms, suggesting that there can be potential issues with clone IDs across the array platforms⁶ ⁷. The most reliable genes selected for validation in this study were those that had the same fold change across the 3 platforms. In addition to selecting candidate genes from the microarray analysis for validation, we included Forkhead Box P3 (FOXP3), Perforin (PFN1) and Granzym B (GRZYB)⁸ ⁹ ¹⁰, as these genes were previously reported to be increased in AR in blood.

Project 2—Verification and Validation of AR Specific Genes by Q-PCR

Candidate Gene Selection Criteria for Q-PCR

The aim of the second project of the study was to verify the results obtained from the first project using Q-PCR, and identify a smaller set of genes that could be used for further validation studies for AR specific blood biomarkers. We selected 32 genes from the 525 overlapping genes based on the following criteria: i) 12 genes were chosen as genes of biological relevance and differential expression across all 3 platforms (up-regulated genes: IFNGR1, RNF130, PBEF1, ITGAX, RYBP, CFLAR, PSEN1; down-regulated genes: NKTR, GOLGA8A, IL32, MCM7 and NFATC3), ii) 17 genes were selected with biological relevance and differentially expression on at least 2 platforms with fold change>1.5 and q<3% (up-regulated genes: DUSP1, SLP1, IL8RAP, STAT3, F2RL1, FCGR1A, IL6R, PTPRC, STAT1, IL1RAP, TLR8, TNFAIP6; down-regulated genes: MAPK9, PHLDA1, PTPN11, ZP70, PLCG1), iii) 3 genes were selected based on a priori published data on their potential as blood based AR biomarkers; perforin or PFN1¹¹, granzyme B or GRZYB¹⁰ ¹¹ and the forkhead box P3 gene, FOXP3⁸.

Verification by Q-PCR

Expression of the 32 genes in AR was verified using TagMan Gene Expression Assay products (Applied Biosystems, Foster City, Calif.) in 34 samples (Verification Set; 17 AR, 17 STA). These samples had been previously run on one or more arrays, and samples were selected based only on the criteria that the remaining sample had adequate RNA quantity and quality. The primer list for each of the 10 genes, and their corresponding ABI probe IDs are provided in FIG. 6A. Expression of the 10 most significant genes are shown in FIG. 6B. For each gene, no fluorescent signal was generated by these assays when genomic DNA was used as a template, which confirms that the assay measured only mRNA.

Ribosomal RNA 18s was used for endogenous RNA and cDNA quality control. 18s was chosen as a control gene, on the basis of this gene having relatively constant gene expression across AR samples and between AR and STA samples on array experiments. In order to control for possible variations among PCRs performed on different days, the expression of all genes were also assessed in Universal Human Reference RNA for each Q-PCR run. The assays were highly reproducible, with coefficient of variation less than 0.18 among 6 runs for each gene assessed in the Universal Human Reference RNA. Q-PCR was carried out in a 384-well plate, with each gene expression measurement performed in duplicates. The relative amount of RNA expression was calculated using a comparative C_(T) method.

In addition to verifying the microarray results, another aim of project 2 was to refine the selection of a smaller, stringent set of genes that could be used validated as biomarkers for AR in independent patient samples. We selected the 10 most significant genes from the 32 genes verified by Q-PCR. Student's t-test was used for establishing the significance of gene expression differences between AR and STA samples by Q-PCR. All differences were considered statistically significant for p<0.05. These 10 genes were next validated for their significant expression differences in AR, by Q-PCR, in 42 independent samples (Validation Set 1; 21 AR, 21 STA) (FIG. 8).

Classification Modeling for AR

To ascertain the best classification model for a minimum gene set for AR diagnosis in blood, we applied 10 different classification models to the 10-gene set¹² he classification models used included Classification ViaRegression, variants of Bayesian models (NaiveBayes, MultinomialNaiveBayes, BayesNet, etc.), multilayer perceptron, logistic regression (multinomial and linear) and sequential minimal optimization. Each of these models were trained using the Q-PCR data for Validation set 1, and re-validated using Q-PCR data from 64 independent samples (Validation Set 2; 32 AR, 32 STA) (FIG. 9). Out of the 10 classification models, multinomial logistic regression model using 5 genes (DUSP1, MAPK9, NKTR, PBEF1, PSEN1) performed the best (sensitivity=100%, specificity=93.8%, PPV=94.1% 1, NPV=100%). FIG. 10 provides comparisons of all classification models used.

A first logistic regression model using the 5-genes set is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {10.0395 + {({1.8364*{DUSP}\; 1})} + {({{- 1.9182}*{MAPK}\; 9})} + {({{- 1.8939}*{NKTR}})} +} \\ {{({4.9901*{PBEF}\; 1})} + {({2.124*{PSEN}\; 1})}} \end{matrix}}}{1 + ^{\begin{matrix} {10.0395 + {({1.8364*{DUSP}\; 1})} + {({{- 1.9182}*{MAPK}\; 9})} + {({{- 1.8939}*{NKTR}})} +} \\ {{({4.9901*{PBEF}\; 1})} + {({2.124*{PSEN}\; 1})}} \end{matrix}}}$

A second logistic regression model using a minimal 5-gene set is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {{- 18.6137} + {({3.8327*{DUSP}\; 1})} + {({{- 5.1524}*{MAPK}\; 9})} + {({{- 3.6056}*{NKTR}})} +} \\ {{({8.3652*{PBEF}\; 1})} + {({3.353*{PSEN}\; 1})}} \end{matrix}}}{1 + ^{\begin{matrix} {{- 18.6137} + {({3.8327*{DUSP}\; 1})} + {({{- 5.1524}*{MAPK}\; 9})} + {({{- 3.6056}*{NKTR}})} +} \\ {{({8.3652*{PBEF}\; 1})} + {({3.353*{PSEN}\; 1})}} \end{matrix}}}$

The logistic regression model using a minimal 5-gene set based on the Stanford samples is shown in the following equation:

$\theta = \frac{^{\begin{matrix} {{- 17.2226} + {({1.0735*{IFNGR}\; 1})} + {({{- 7.7795}*{MAPK}\; 9})} + {({4.2314*{PBEF}\; 1})} +} \\ {{({3.7546*{PSEN}\; 1})} + {({2.9103*{RNF}\; 130})}} \end{matrix}}}{1 + ^{\begin{matrix} {{- 17.2226} + {({1.0735*{IFNGR}\; 1})} + {({{- 7.7795}*{MAPK}\; 9})} + {({4.2314*{PBEF}\; 1})} +} \\ {{({3.7546*{PSEN}\; 1})} + {({2.9103*{RNF}\; 130})}} \end{matrix}}}$

In each equation above, θ is the predicted probability. Based on the Receiver Operating Characteristic (ROC) curve, a cutoff of θ=0.5 was selected to determine whether the predictive class was an AR or STA phenotype.

Querying the Influence of Demographic and Clinical Confounders on the 5 Gene-Set

Data was collected on 18 demographic and clinical variables on all samples processed for Q-PCR validation studies. Data was provided for SNSO1 samples by PPD (Brandon Graham, Jason Berg, David Ikle). The following 18 variables were studied: post-transplant time, recipient age, recipient gender, donor gender, donor source, donor age, cytomegalovirus/CMV, Epstein-Barr virus/EBV and BK virus/BKV viremia at sampling, concomitant bacterial infection, presence of donor specific antibodies (DSA) or panel reactive antibodies (PRA), usage of induction therapy, maintenance steroid usage, use of different calcineurin inhibitor drugs (FK506 or cyclosporine A/CSA), use of different anti-metabolites (mycophenolate mofetil/MMF or azathioprine/AZA), the white blood cell count (WBC) and the hematocrit (HCT) levels at the time of blood sampling (FIG. 11). Pearson correlation coefficients were calculated to determine any influence of 18 different clinical parameters as potential confounders for expression of the 5-genes set (FIG. 12). This analysis showed that the expressions of the 5 gene-set for diagnosis of AR, was not affected by any of the 18 confounders examined (maximum |r|=0.27). Additionally, t-test was performed for 13 clinical confounders with categorical values (recipient gender, donor gender, donor source, steroid-free versus steroid-based, CMV viremia, EBV viremia, BK viremia, bacterial infection, DSA, PRA, use of induction therapy or not, use of FK506 versus CSA, use of MMF versus AZA).

Querying the Influence of AR Severity (BANFF Scores), AR Type (Humoral or Cellular Rejection), Peritubular C4d Positivity at AR, on the Expression of the 5 Gene-Set

Next, any associations between Banff scores, peri-tubular C4d+/− at biopsy and the presence of either humoral or cellular rejection, and expression of the 5-gene set for AR diagnosis, were examined. As most samples had pared biopsies for evaluation, samples were divided into 2 groups according to Banff score grades (1=Banff scores≦1A; 2=Banff score≧1B); C4d⁺ was determined if there was ≧10% of peritubular capillaries (PTC) had C4d deposits. No correlation was found between AR probability prediction by the 5 gene-set and Banff scores (r=−0.13), humoral versus cellular rejection (r=0.07), and presence or absence of peri-tubular C4d+ (r=0.2) (FIG. 14).

Project 3—Prediction of AR Injury by the 5 Gene-Set, Prior to Indicated Graft Biopsy and Graft Dysfunction

The aim of the third project in this study was to test whether the model with 5-gene set can predict the onset of AR, prior to any graft dysfunction or indicated biopsy. We applied our regression model to samples sequentially collected 1 to 6 month prior and after a biopsy proven AR episode. The samples were grouped into the following categories: i) samples collected 3 to 6 months prior to the onset of AR, ii) samples collected 0 to 3 months prior to AR, iii) samples collected at the time of AR, iv) samples collected 0 to 3 month after the AR episode, v) samples collected 3 to 6 months after the AR episode. FIG. 15 provides mean AR probabilities along with standard errors and calculated creatinine clearances¹³ for each group of samples.

Biological Relevance Analysis of the 5 Gene-Set

To examine the key biological processes that could influence the peripheral transcriptome in AR, we applied the 525 cross platform AR specific gene-set for biological relevance analysis using the PANTHER classification system (http(colon)//www(dot)pantherdb(dot)org). These molecular pathways suggest active protein trafficking/transport and lipid metabolism in key pathways involved in immunity, cell cycling, apoptosis and differentiation. Ingenuity Pathway Analysis (http(colon)//www(dot)ingenuity(dot)com) shows that these genes drive canonical pathways in involved in apoptosis, and regulated by the JAK/STAT pathway (FIG. 16).

REFERENCES FOR SUPPLEMENTAL METHODS

1. Racusen L C, Solez K, Colvin R B, et al. The Banff 97 working classification of renal allograft pathology. Kidney Int 1999; 55(2):713-23.

2. Li L, Ying L, Naesens M, et al. Interference of globin genes with biomarker discovery for allograft rejection in peripheral blood samples. Physiol Genomics 2008; 32(2):190-7.

3. Sarwal M, Chua M S, Kambham N, et al. Molecular heterogeneity in acute renal allograft rejection identified by DNA microarray profiling. N Engl J Med 2003; 349(2):125-38.

4. Efron B, Tibshirani R. Empirical bayes methods and false discovery rates for microarrays. Genet Epidemiol 2002; 23(1):70-86.

5. Chen R, Li L, Butte A. AILUN: reannotating gene expression data automatically. Nature Methods 2007; 4(11):In Press.

6. Perez-Iratxeta C, Andrade M A. Inconsistencies over time in 5% of NetAffx probe-to-gene annotations. BMC Bioinformatics 2005; 6:183.

7. Harbig J, Sprinkle R, Enkemann S A. A sequence-based identification of the genes detected by probesets on the Affymetrix U133 plus 2.0 array. Nucleic Acids Res 2005; 33(3):e31.

8. Li B, Hartono C, Ding R, et al. Noninvasive diagnosis of renal-allograft rejection by measurement of messenger RNA for perforin and granzyme B in urine. N Engl J Med 2001; 344(13):947-54.

9. Dugre F J, Gaudreau S, Belles-Isles M, Houde I, Roy R. Cytokine and cytotoxic molecule gene expression determined in peripheral blood mononuclear cells in the diagnosis of acute renal rejection. Transplantation 2000; 70(7):1074-80.

10. Simon T, Opelz G, Wiesel M, Ott R C, Susal C. Serial peripheral blood perforin and granzyme B gene expression measurements for prediction of acute rejection in kidney graft recipients. Am J Transplant 2003; 3(9):1121-7.

11. Vasconcellos L M, Schachter A D, Zheng X X, et al. Cytotoxic lymphocyte gene expression in peripheral blood leukocytes correlates with rejecting renal allografts. Transplantation 1998; 66(5):562-6.

12. Witten I H, Frank E. Data Mining: Practical machine learning tools and techniques. 2nd ed. San Francisco: Morgan Kaufmann; 2005.

13. Schwartz G J, Haycock G B, Edelmann C M, Jr., Spitzer A. A simple estimate of glomerular filtration rate in children derived from body length and plasma creatinine. Pediatrics 1976; 58(2):259-63.

Example II Methods

A total of 40 time and immunosuppression matched peripheral blood samples from liver and heart transplant patients were selected for 384-well format multiplex ABI-TaqMan qPCR validation across the 10 gene-set identified in Example I. 16 of the samples were from liver transplant patients with 8AR and 8 stable without rejection (STA) and 24 samples were from heart transplant patients (9AR; 10STA); for virus confounder control, 5 STA samples with Cytomegalovirus infection (CMV) were included. The significance differential gene expression was determined by T test, with p<0.05 considered as significant.

Results

For the same 10 gene-set validation in liver transplant, 2 of the 10 genes (PSEN1 and RYBP) were highly significant for confirming peripheral blood AR diagnosis (p=0.05; p=0.04) in comparison to STA; the remaining 8 genes showed the same-direction fold changes as seen in kidney transplant, with no statistically significant difference between AR and STA. A prediction model generated from the 10 gene-set derived from peripheral blood of kidney transplant (ROC=97.6%) was applied to liver sample PCR dataset for AR classification, and revealed a prediction score of 86% sensitivity, 75% specificity, 88% PPV and 78% NPV, with 2 AR and 1 STA of the 16 samples misclassified. For 10 gene-set validation in heart samples, 3 of the 10 genes (CFLAR, ITGAX and NKTR) were differentially expressed in AR compared to STA (p=0.04; p=0.05; p=0.01). Interestingly, 1 of these 3 genes (CFLAR) shows significant opposite-direction fold changes (p=0.04, fold change=−2.2) as seen in kidney AR samples previously (p=0.002, fold change=3.1). Furthermore, 2 of the 3 genes were significantly down-regulated in CMV when compared to AR. There was no overlap in the significant differentially expressed genes between liver and heart. The identity of the 10 genes tested is shown in the previous discussion with kidney transplant samples.

Conclusion

A highly specific and sensitive biomarker of 10 genes derived and tested by cross-platform microarray analysis in kidney transplant can be used to predict AR in liver transplant. 2 of the 10 genes were differentially expressed in AR for both liver and heart, but no overlap of significantly differentially regulated genes between liver and heart suggests that rejection injury not only shares common biological pathways among different organ transplants, but also exhibits tissue specificity.

The 10 genes tested can be used to diagnose AR in both liver and heart transplant rejection, regardless of rejection grade.

Although the foregoing invention has been described in some detail by way of illustration and example for purposes of clarity of understanding, it is readily apparent to those of ordinary skill in the art in light of the teachings of this invention that certain changes and modifications may be made thereto without departing from the spirit or scope of the appended claims.

Accordingly, the preceding merely illustrates the principles of the invention. It will be appreciated that those skilled in the art will be able to devise various arrangements which, although not explicitly described or shown herein, embody the principles of the invention and are included within its spirit and scope. Furthermore, all examples and conditional language recited herein are principally intended to aid the reader in understanding the principles of the invention and the concepts contributed by the inventors to furthering the art, and are to be construed as being without limitation to such specifically recited examples and conditions. Moreover, all statements herein reciting principles, aspects, and embodiments of the invention as well as specific examples thereof, are intended to encompass both structural and functional equivalents thereof. Additionally, it is intended that such equivalents include both currently known equivalents and equivalents developed in the future, i.e., any elements developed that perform the same function, regardless of structure. The scope of the present invention, therefore, is not intended to be limited to the exemplary embodiments shown and described herein. Rather, the scope and spirit of present invention is embodied by the appended claims. 

What is claimed is:
 1. A method of monitoring a subject who has received an allograft for an acute rejection (AR) response, said method comprising: (a) evaluating the expression level of at least one gene in a sample from said subject to obtain a gene expression result, wherein said at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 and PSEN1; and (b) employing said gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in said subject, thereby monitoring said subject for an AR response.
 2. The method according to claim 1, wherein the expression level of NKTR, MAPK9, DUSP1, PBEF1 and PSEN1 is assessed.
 3. The method according to claim 1, wherein said evaluating step further comprises evaluating the expression level of one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP.
 4. The method according to claim 1, wherein said evaluating step comprises assaying said sample for an expression product of said at least one gene.
 5. The method according to claim 4, wherein said expression product is a nucleic acid transcript.
 6. The method of claim 5, wherein said assessing step comprises an RT-PCR assay.
 7. The method of claim 6, wherein said RT-PCR assay is a quantitative RT-PCR assay.
 8. The method according to claim 1, wherein said sample is a blood sample.
 9. The method according to claim 1, wherein said allograft is a kidney allograft.
 10. The method according to claim 1, wherein said gene expression result is employed to predict the occurrence of an AR response between 6 and 3 months in advance.
 11. The method according to claim 1, wherein said method further comprises evaluating the expression level of at least one gene selected from Tables 3 or 4 in said sample from said subject to obtain a second gene expression result and employing said second gene expression result to determine whether said subject has a graft tolerant phenotype.
 12. A system for monitoring a subject who has received an allograft for an acute rejection (AR) response, said system comprising: (a) a gene expression evaluation element configured for evaluating the expression level of at least one gene in a sample from said subject to obtain a gene expression result, wherein said at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 and PSEN1; and (b) a phenotype determination element for employing said gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in said subject, thereby monitoring said subject for an AR response.
 13. The system according to claim 12, wherein said gene expression evaluation element comprises at least one reagent for assaying a sample for an expression product of said at least one gene.
 14. The system according to claim 13, wherein said expression product of said at least one gene is a nucleic acid transcript.
 15. The system according to claim 14, wherein said gene expression evaluation element comprises PCR primers specific for said at least one gene.
 16. The system according to claim 12, wherein the expression level all of NKTR, MAPK9, DUSP1, PBEF1 and PSEN1 is assessed.
 17. The system according to claim 12, wherein said gene expression evaluation element is further configured for evaluating the expression level of one or more of: CFLAR, RNF-130, IFNGR1, ITGAX and RYBP.
 18. The system according to claim 12, wherein said phenotype determination element comprises a reference expression value for said at least one gene.
 19. A kit for monitoring a subject who has received an allograft for an acute rejection (AR) response, said kit comprising: (a) a gene expression evaluation element for evaluating expression of at least one gene in a sample to obtain a gene expression result, wherein said at least one gene is selected from the group consisting of: NKTR, MAPK9, DUSP1, PBEF1 and PSEN1; and (b) a phenotype determination element for employing said gene expression result to predict the onset of an AR response, to diagnose an AR response, and/or to characterize an AR response in said subject, thereby monitoring said subject for an AR response.
 20. A computer program product for monitoring a subject who has received an allograft for an acute rejection (AR) response, wherein said computer program product, when loaded onto a computer, is configured to employ a gene expression result from a sample derived from said subject to determine an AR monitoring result and provide said AR monitoring result to a user in a user-readable format, wherein said AR monitoring result is selected from predicting the onset of an AR response, diagnosing an AR response, and/or characterizing an AR response in said subject, wherein said gene expression result comprises expression data selected from one or more of: NKTR, MAPK9, DUSP1, PBEF1 and PSEN1. 